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Abstract 

There is good evidence that full general relativity is non-integrable or even chaotic. We point 
out the severe repercussions: differentiable Dirac observables and a reduced phase space do not exist 
in non-integrable constrained systems and are thus unlikely to occur in a generic general relativis¬ 
tic context. Instead, gauge invariant quantities generally become discontinuous, thus not admitting 
Poisson-algebraic structures and posing serious challenges to a quantization. Non-integrability also 
renders the paradigm of relational dynamics cumbersome, thereby straining common interpretations 
of the dynamics. We illustrate these conceptual and technical challenges with simple toy models. 
In particular, we exhibit reparametrization invariant models which fail to be integrable and, as a 
consequence, can either not be quantized with standard methods or lead to sick quantum theories 
without a semiclassical limit. These troubles are qualitatively distinct from semiclassical subtleties in 
unconstrained quantum chaos and can be directly traced back to the scarcity of Dirac observables. As 
a possible resolution, we propose to change the method of quantization by refining the configuration 
space topology until the generalized observables become continuous in the new topology and can 
acquire a quantum representation. This leads to the polymer quantization method underlying loop 
quantum cosmology and gravity. Remarkably, the polymer quantum theory circumvents the prob¬ 
lems of the quantization with smooth topology, indicating that non-integrability and chaos, while a 
challenge, may not be a fundamental obstruction for quantum gravity. 


1 Introduction 

In order to extract physical information from gauge theories, one either constructs gauge invariant ob¬ 
servables or fixes a gauge. In a canonical language, this usually means either constructing functions which 
Poisson-commute on the constraint surface with the gauge generating first class constraints or to fix the 
flow of the latter by singling out points within each gauge orbit l||2]. Such gauge invariant observables 
are thus the ‘constants of motion’ of the gauge generators and referred to as Dirac observables. In the 
quantum theory, one typically deals with gauge symmetries by canonically quantizing according to the 
Dirac (and sometimes even the reduction) method or by performing a Faddeev-Popov type gauge fixing 
in the corresponding path integral formulation. In either method, Dirac observables assume a pivotal 
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role in the interpretation of the quantum theory. The reason for this is that the operators on the physi¬ 
cal Hilbert space, obtained after solving the (quantum) constraints, are the quantum Dirac observables. 
Also, if one uses a path integral quantization, the path integral acts as a projector onto the physical 
Hilbert space [3,[4|, and (quantum) Dirac observables are again essential to characterize the resulting 
quantization. 

Among the gauge theories, general relativity possesses special attributes: it features spacetime diffeo- 
morphism symmetry and, in its canonical formulation, is a totally constrained system j5]-[7], which means 
that also its dynamics is generated by a (first class) Hamiltonian constraint. This has deep consequences: 
Firstly, diffeomorphism invariant observables are necessarily non-local in a spacetime manifold sense (i.e., 
they must involve spatial derivatives of infinitely high order) [8] which poses challenges to connecting them 
to the (background dependent) local observables of quantum field theory 9,10 . Secondly, gravitational 
Dirac observables are coordinate and, in particular, time coordinate independent; they are ‘constants of 
motion’ of the Hamiltonian constraint, thus appearing to be non-dynamical, which leads to the infamous 
‘problem of time’ in the quantum theory 


11 13 


A possible resolution to the seemingly non-local and non-dynamical character of gravitational observ¬ 
ables is often hoped to be found within the relational nature of general relativity itself: localization and 
evolution of systems and dynamical quantities are only meaningful relative to other dynamical quantities 
- rather than with respect to a background. Gravitational observables capturing the relative localiza¬ 


tions and dynamics should thus be relational objects that correlate different dynamical quantities 14 ■ 20 


Such relational Dirac observables are often also called ‘evolving constants of motion’. Classically, much 
progress has been made in defining and, at least formally, constructing relational observables 20-26 


In particular, there now exist perturbative expansions, applicable to fixed backgrounds and to pertur¬ 
bations around symmetry reduced sectors (e.g. in cosmology), for constructing approximate observables, 
which in their lowest order approximation reproduce the local observables of field theory on a fixed back¬ 
ground 27 28 . On the other hand, much work remains to be done in the quantum theory - for both 


the finite dimensional and the field theory cases - where no general method for representing the quan¬ 
tum analogues of relational observables exists. As a possible remedy, effective techniques for evaluating 
and constructing relational observables in the semiclassical limit of finite dimensional systems have been 
developed which, in principle, can also be generalized beyond the semiclassical regime and to the field 


theory case 29-32 


Unfortunately, gravitational Dirac observables (relational or otherwise) are notoriously hard to find 
explicitly as their construction requires solving the dynamics (see [20||26| for recent reviews). In fact, 
for full general relativity no Dirac observables are known with the exception of the ADM charges in 
asymptotically flat spacetimes 6 , the boundary charges in asymptotically anti de Sitter spacetimes 


33-35 or, perhaps, the deparametrized relational observables in dust filled spacetimes 136,37 . However, 


these exceptions do not account for a generic situation in full general relativity since the first two examples 
require special asymptotic symmetries, while the third relies on a specially engineered matter content. 
Indeed, Anderson and Torre showed that, except for the asymptotic symmetries, there do not exist any 


further (generalized) symmetries in full general relativity 38 


Relational Dirac observables, in particular, are difficult to construct and evaluate explicitly because 
they require internal time or ‘clock’ variables which parametrize the orbits of the gauge generating 
constraints. Since fixing the ‘clocks’ to a particular set of values should single out a point on each gauge 
orbit, the existence of such ‘clock’ variables is directly connected to the existence of good gauge fixing 
conditions for the flows of the constraints. Indeed, relational observables can be viewed as gauge invariant 
extensions of gauge restricted quantities ||2 21,22 . In complete analogy to the Gribov problem in non- 
abelian gauge theories, however, there will generically exist global obstructions to the existence of ideal 
‘clocks’ whose level surfaces are intersected precisely once by every classical trajectory; this has become 
known as the ‘global problem of time’ 11 12 16 23 24 29 -32 39 41 . Furthermore, a generic ‘clock’ 


variable will couple to the ‘evolving’ degrees of freedom which severely complicates the resolution of the 
quantum relational dynamics 9,3l|42 43 . In other words, ideal ‘clocks’ will in general not exist. This is 
emphasized in the work by Torre [44[ ) showing that general relativity is not a parametrized system, i.e. 
there does not exist an ideal (set of) clock variables. This is an important point since - as we shall see - 
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the problem we are discussing here does not arise for parametrized systems. 

In (loop) quantum cosmology, one usually turns a blind eye to this challenge by employing highly 
idealized and decoupled matter ‘clocks’ such as a free, massless scalar 45 46 or dust field 47 , or by 


resorting to model specific geometric internal time variables 48 . While this is useful for gaining insights 


into non-generic quantum cosmological models, it clearly does not capture the qualitative features of 
relational dynamics and observables in a generic cosmological model. The effective techniques 29-32 


for evaluating quantum relational dynamics have been developed precisely to remedy the restriction to 
special ‘clocks’ and, instead, to generalize to arbitrary degrees of freedom as internal times. In particular, 
this method permits one (a) to switch the ‘clock’ in the quantum theory, thereby enabling one to evolve 
past ‘clock’ pathologies by patching together evolutions in different internal times, and (b) to construct 
and give meaning to transient/local observables, both of which seem indispensable in the light of the 
‘global problem of time’. 

However, as argued in 31 ^ the ‘global problem of time’ will become unsurmountable in the pres¬ 
ence of chaos. Chaos generically prevents the existence of differentiable internal time variables which 
parametrize the orbits of the first class constraints (even locally) in every open neighbourhood of the 
constraint surface. (It does not prevent the existence of discontinuous internal time variables which, 
however, cannot be inserted in Poisson brackets or the existence of local ‘clocks’ which parametrize the 
orbits away from chaotic regions.) This has severe repercussions for quantum relational dynamics which, 
even semiclassically, seems to break down generically 31 . The absence of a good relational clock also 


implies that the Gribov problem becomes generically untractable (at least in some neighbourhoods of the 
constraint surface) in the presence of chaos; a given gauge might either miss ‘most’ trajectories and/or 
have densely many solutions (see also the main body of this article). Since chaos features in generic dy¬ 
namical systems [49-151 , the understanding of relational dynamics and observables, gained from simple, 


solvable cosmological models, will thus, in fact, be misleading for a general cosmological model, let alone 
full quantum gravity. Under such circumstances, it seems advantageous to consider physical transition 
amplitudes (involving the projector onto the constraints) between kinematical states 19 52 , instead of 


the more cumbersome relational dynamics. Indeed, such transition amplitudes are more general and, in 
principle, admit a useful dynamical interpretation also when relational dynamics breaks down. 

This article is devoted to better understanding the repercussions of chaos — or, more generally, non- 
integrability - for the existence of Dirac observables, for relational dynamics and the quantization of 
constrained systems. We shall argue that, as a consequence of non-integrability which implies the ab¬ 
sence of sufficiently many smooth ‘constants of motion’, it is not merely a technical problem which 
inhibits the construction of Dirac observables for full general relativity or even a generic cosmological 
model, but a problem in principle: Dirac observables for full general relativity (without special asymptotic 
symmetries or matter content) almost certainly do not exist. Before expanding on this in the main body 
and until then causing an unnecessary confusion, we emphasize that this does not imply that gauge 
(i.e., diffeomorphism) invariant quantities generically do not exist. Quite the contrary: gauge invariant 
functions on the constraint surface which are constant on the gauge orbits will always exist0 However, 
these will generically be discontinuous (or singular, see 24 25 53 and main body of this article), i.e., 


non-differentiable and thus not Dirac observables in the literal sense since they cannot be inserted into 
Poisson brackets with other observables. That is, these generalized gauge invariant quantities will not 
admit any Poisson algebraic structure. We emphasize that with the term ‘Dirac observable’ we refer to 
phase space functions which are smooth (or sufficiently differentiable in all directions) on the constraint 
surface. 

What might appear as an innocent observation, has deep consequences, in particular for the quantum 
theory: Firstly, it implies that a reduced phase space does not exist under such circumstances since the 
Poisson (or rather Dirac) bracket between any two functions on the reduced phase space - if it existed 
- are equivalent to a Poisson bracket between two Dirac observables on the constraint surface [2j. But 
such differentiable Dirac observables are likely not to exist generically. This suggests that a reduced phase 
space for full general relativity (without special asymptotic symmetries or matter content) almost certainly 


1 See especially Sec. IV C 4 of [31 . 

2 Such functions can, for instance, be defined by averaging non-invariant functions over gauge orbits., 
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does not exist either. 

Secondly, in the absence of Dirac observables one lacks an algebraic structure which can be promoted 
into the quantum theory. This directly rules out reduced quantization. While the failure of the reduction 
method by itself is not a great surprise, the absence of Dirac observables and algebraic relations also 
poses a challenge to a Dirac or path integral quantization through standard methods (if possible at all) 
because it will become very non-trivial to relate the observables on a physical Hilbert space to those of 
the classical theory. This severely complicates a semiclassical limit and the interpretation of the quantum 
theory. 

Thirdly, as we shall illustrate by means of toy models, the absence of sufficiently many Dirac observ¬ 
ables indeed seems to generically prevent the very existence of a semiclassical limit - both in terms of 
physical states and transition amplitudes - of a quantum theory obtained by a standard Dirac procedure. 
This supports earlier reports in the literature of an absence of good semiclassical states in a chaotic 
cosmological model 31 54 . We emphasize at this point already that these troubles are qualitatively 


and conceptually distinct from the usual semiclassical subtleties in unconstrained quantum chaos where 
typically at least classically periodic orbits can be reproduced and wave packets have at least a short 
term coherence 49-51 . In the present constrained case, the absence of semiclassical states is directly 


tied to the absence of sufficiently many Dirac observables. Such Dirac quantized models can generically 
not reproduce any of the classical dynamics and may thus be considered iii0 In particular, we shall 
also exhibit a system which has meaningful classical dynamics, yet can neither be quantized through 
the reduced nor the standard Dirac method, thereby providing an extreme example for a sick quantum 
theory (and, in particular, the absence of a semiclassical limit). 

All these observations seem, at first sight, to be rather devastating news for the quantization of non- 
integrable constrained systems in general and general relativity in particular. On a second look, however, 
these properties turn out to be less dramatic and do not necessarily constitute a fundamental road 
block to quantum gravity. Indeed, one has to remember that quantizing a theory is a highly ambiguous 
procedure. These findings simply force us to change our preferred method - and, more precisely, topology 
- for quantization. The mentioned troubles with the standard techniques can be directly traced back to 
the discontinuity of the observables which impedes their quantum representations. However, the notion 
of (dis)continuity requires a topology and the generalized observables are discontinuous in the standard 
(Euclidean metric) topology. The key idea we shall propose here for resolving these troubles and usefully 
translating these observables into a quantum theory will be to refine the topology until the generalized 
observables become continuous in the new topology. 

This leads to the good news of this article: Introducing such a refined topology leads to an interesting 
quantum theory with sufficiently many observables in all presented examples where the standard methods 
produce a failure. 

Clearly, there is a price to pay for curing the illnesses of standard quantization methods: (a) One has 
to deal with an unfamiliar discrete topology and the ensuing technical consequences, (b) one obtains a 
novel parameter in the quantum theory measuring the step size of translation operators which requires a 
physical interpretation, and (c) the spectra of the conjugate observables will become (Bohr) compactified. 
However, there already exists a well-developed machinery for dealing with this alternative quantization 
method. Namely, interestingly, this method is precisely the quantization method of loop quantum cos¬ 
mology (45 46 55 , sometimes known as ‘polymer quantization’ 56-58 , and also appears in a recent 


novel representation of loop quantum gravity 59 . In particular, in loop quantum cosmology, the novel 


parameter (b) introduced through quantization is closely related to the Planck length and thus carries 
an immediate physical interpretation. The results presented in this manuscript are a remarkable success 
of polymer quantization and provide evidence that this method might generally resolve the challenges 
posed by non-integrability, although the task of finding a useful interpretation of the quantum dynamics 
(relational or in terms of transition amplitudes) remains. 

Given the importance of Dirac observables and relational dynamics for constrained systems and their 
quantization in general - and general relativity in particular - and given the deep consequences of non- 


3 This is also distinct from systems such as a hydrogen atom where states of low energy do not reproduce the (wrong) 
classical low energy dynamics, yet highly energetic states come arbitrarily close to the classical high energy case. 
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integrability featuring in generic systems, it is surprising that this issue has been largely ignored in the 

which, however, only briefly touch 


recent literature. The exception are older works 24 25 38 53 60-62 


on the existence of classical Dirac observable^ and are rather brief on the consequences for the quantum 
theory^] The first concrete analysis of the deep ramifications of chaos for the quantum theory of con¬ 
strained systems and relational dynamics in particular has only been carried out recently 31 . However, 
in our opinion, the consequences of non-integrability for quantum gravity cannot be underestimated and 
warrant further study. In this spirit, we hope that the present work opens up a new avenue of research 
on chaos and non-integrability in (quantum) gravity]^ 

The remainder of this manuscript is organized as follows. In section [2j we shall make the notion 
of non-integrability more precise and discuss the evidence for non-integrability in general relativity. In 
order to display a concrete analysis of the consequences of non-integrability, we introduce a toy model in 
section [3] with ergodic dynamics in configuration space. This model, while not being fully non-integrable 
or chaotic, lacks sufficiently many Dirac observables and is shown not to possess a reduced phase space. 
The advantage of this simple ‘semi-integrable’ model is that it can still be treated analytically (in contrast 
to most chaotic models), yet features generic behaviour of non-integrable models. We quantize this model 
by means of the Dirac method and standard topology in section [3J exhibiting its failure to produce a 
good semiclassical limit. Thereupon, in section [5] we refine the topology and polymer quantize the same 
model, showing how an interesting quantum theory with sufficiently many observables arises. In section 
[6j we discuss three further twists on this toy model, illustrating various degrees of integrability (however, 


none of which leads to full non-integrability). In particular, in section 6.1 we discuss a model whose 


classical trajectories are all ergodic and which is shown to prevent a quantization by standard methods, 
yet which permits a polymer quantization. We conclude with a discussion in section [7] 


2 Non-integrability and Dirac observables 

It is useful to review the concepts of integrability and non-integrability in unconstrained Hamiltonian 
systems. We shall do so for finite dimensional autonomous systems. 

Consider a system described by a 2iV-dimensional phase space V and time independent Hamiltonian 
H. The system is called integrable if there exist N independent]^] constants of motion T — {Fj,..., Fy} 
on V with H £ J- which are in involution 


{F k ,Fi}= 0 


( 2 . 1 ) 


and thus, by Frobenius’ theorem, integrate to a V-dimensional submanifolds Mp of V, determined by 
the values of F\,... ,Fjy. Under mild conditions, one can show that Mp ~ T fe x k where T fc is a 
fc-dimensional torus 64 . That is, in an integrable system, the classical trajectories lie on IV-dimensional 


submanifolds Mp of phase space. 

By contrast, an unconstrained Hamiltonian system is called non-integrable if no smooth (or differen¬ 
tiable) constants of motion exist globally on V which are independent of H. In this case, trajectories are 
not confined to an iV-dimensional submanifold of V, but only to a (2 N — l)-dimensional energy surface; 
a typical trajectory will explore a region of co-dimension 1 in phase space. Non-integrability features 
in generic unconstrained Hamiltonian systems 


49-51 


conditions for the absence of smooth constants of motion 64 


and there exist concrete theorems stating precise 
There are various characterizations of 


non-integrable systems 49-51 64 ; e.g., they may be ergodic or chaotic. Ergodic or chaotic systems are 


non-integrable, however, the converse is not necessarily true. The precise definition of chaos or ergodicity 
are unimportant for our purposes. What matters is the more general notion of non-integrability and the 
associated absence of constants of motion. 


4 In the older literature, Dirac observables are sometimes also referred to as ‘perennials’. 

5 In particular [53 notes that (smooth) Dirac observables will not exist for an ergodic dynamics. Our examples will also 
be ergodic at lea st in configuration space, resulting in the lack of sufficiently many Dirac observables. 

6 See also |63| for an interesting application of chaos and random matrix theory in quantum gravity to provide evidence 
that the spectrum of the volume operator features a gap at the zero eigenvalue. 

' That is, the differentials of the F-i are linearly independent everywhere on phase space. 
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In analogy to these notions for unconstrained Hamiltonian systems, we shall now define three distinct 
notions of (non-)integrability for constrained systems. Again, we restrict to the finite dimensional case; 
the field theory case might be defined analogously in terms of local degrees of freedom. Consider a 
system described by a 2IV-dimensional kinematical phase space V and m\ irreducible first class and 2 m 2 
irreducible second class constraints such that - for a solvable system - one would expect 2(N — mi — m 2 ) 
independent physical degrees of freedom. Denote by C the constraint surface and by 7 a (copy of the) 
gauge orbit. 

We begin with a physically somewhat superfluous notion of (non-)integrability, however, mention it 
for technical completeness. 


strong (non-)integrability: A constrained system shall be referred to as strongly integrable if there 
exist N independent functions Fi,..., Fn on V, containing the constraints, which are in involution 
(2.1) in terms of strong equalities. These functions, again, integrate to iV-dimensional hypersurfaces 
in the kinematical phase space. 


A constrained system shall be called strongly non-integrable if there do not exist any functions, 
independent of the constraints, which strongly commute with the constraints and if the constraints 
themselves do not commute strongly. 


Clearly, a system with second class constraints can never be strongly integrable. Regardless, strong 
equalities are usually physically not very interesting for constrained systems such that we shall not 
further elaborate on this case. 

Instead, of most interest to us will be the notion of (non-)integrability in terms of the physically 
relevant weak inequalities. This is also the notion pertinent to general relativistic systems. 


weak (non-)integrability: A constrained system will be called weakly integrable if (a) there exist 2 (N— 
mi — m 2 ) independent Dirac observables which are also independent of the constraints, and (b) 
N — mi — m 2 of the Dirac observables are weakly in involution, (a) implies that the system 
is reducible, i.e. (each connected component of) C/ 7 constitutes a reduced phase space Vr with 
symplectic structure (see the comments in the introduction and !)■ Moreover, (b) implies that, 
together with the mi first class constraints, the N — m\ — m 2 Dirac observables in weak involution 
form a set of N — m 2 functions in weak involution. By Frobenius’ theorem, these integrate to 
( 7 V+m 2 )-dimensional submanifolds of C. If the generator of dynamics is first class (this is generically 
the case [ 2 ]), a classical trajectory will remain restricted to such a submanifold. 


On the other hand, a constrained system will be called weakly non-integrable if there exist no Dirac 
observables other than the constraints. As argued in the introduction, this entails that a weakly 
non-integrable system is also necessarily irreducible and does not possess a reduced phase space 
(see also 


24 25 for a related discussion of irreducibility). A Marsden-Weinstein type symplectic 


reduction thus fails in such systems. Furthermore, a classical trajectory will generically not remain 
restricted to a (A r +m 2 )-dimensional submanifold of C, instead exploring regions of larger dimension. 


We emphasize that the weak notion of integrability is indeed mathematically weaker than the strong 
notion as it requires less stringent conditions. 

Finally, there is an obvious notion of (non-)integrability for reducible systems. 

reduced (non-)integrability: A constrained system shall be termed reducibly integrable if the system 
is ( 1 ) weakly integrable and thus reducible, and ( 2 ) if the reduced system is integrable on the 
reduced phase space Vr, i.e. if there exist Nr functions Fj,... ,Fn r in involution with respect to 
the Dirac bracket on the 2IVft-dimensional Vr. 

The constrained system will be referred to as reducibly non-integrable if the system is (1) weakly 
integrable and thus reducible, and (2) if the reduced system is non-integrable on Vr, i.e. there exist 
no smooth constants of motion of the surviving physical Hamiltonian on Vr (which are independent 
of this Hamiltonian). 


6 





For example, consider a non-integrable or even chaotic unconstrained Hamiltonian system subject to 
a (time independent) Hamiltonian H on a 27V-dimensional phase space V. Promoting the time t to a 
dynamical variable and thus extending the phase space to a 2 (N + l)-dimensional phase space produces 
a system governed by the Hamiltonian constraint C = pt + H. This system is clearly reducible because 
t = const constitutes a globally valid gauge fixing condition which gives rise to a global Dirac bracket. 
However, since the reduced phase space is simply the original unextended phase space, Vr = V , the 
system is clearly reducibly non-integrable. (See also [24p5| for a related discussion.) On the other hand, 
if the unconstrained system was integrable, the parametrized analogue of it would evidently be reducibly 
integrable. 

In the sequel, we shall focus almost exclusively on weak (non-)integrability as this implies conse¬ 
quences for the existence of Dirac observables, the existence of a reduced phase space and quantization. 
By contrast, we emphasize that there is no fundamental problem in quantizing reducible systems with 
standard techniques either by the Dirac or reduction method as they possess a reduced phase space and 
thus sufficiently many Dirac observables which generally admit algebraic representations in the quantum 
theory. Clearly, it might be technically difficult to solve the quantum dynamics of reducibly non-integrable 
systems but there is no problem in principle. For instance, for parametrized non-integrable systems the 
physical Hilbert space is isomorphic to the Hilbert space of the original unparametrized system. It will be 
a challenge to exactly solve the Schrodinger (constraint) equation, but there exist useful approximation 


techniques and other powerful methods to extract physical information from such systems 49 - 51 


There is strong evidence that a generic general relativistic system fails to produce Dirac observables 
and thereby constitutes a weakly non-integrable or even chaotic system. (For the field theory case, the 
absence of Dirac observables will be a sufficient definition of weak non-integrability too.) The precise 
definition of chaos for general relativity is unimportant for our purposes. However, we note that there 
exist appropriate definitions of chaos for generally covariant systems in terms of defocusing and ergodicity 
of classical trajectories which does not require the notion of exponential divergence 65 or in terms of the 


topological entropy which measures the complexity of the solution space through the number of closed 
orbits 49 50,66-69 . What typically happens in systems which are chaotic according to such definitions 


is that a generic classical trajectory explores regions of large dimension in the constraint surface - in 
conflict with the notion of weak integrability above. 

We shall list some of the evidence for weak non-integrability of general relativity in increasing strength. 


1 . 


A generic dynamical system is chaotic and this becomes more pronounced the more intricate its 
dynamics 49-51 . Simply by that reasoning - and in light of the complicated general relativistic 


dynamics - it would be a surprise if full general relativity was integrable in any of the above senses. 


2. The N > 3 body problem in Newtonian gravity is chaotic 64 . Given that general relativity 


3. 


reproduces Newtonian gravity in the appropriate limit but is otherwise technically more intricate, 
it is hardly imaginable that general relativity remedies the situation. 

Already one of the simplest cosmological models, namely the closed FRW model with minimally 
coupled massive scalar field, has been shown to be chaotic ![66,70 m (see also |3l] for a discussion 
of relational dynamics in this model). There is analytical and numerical evidence that the space of 
solutions of this model features a fractal structure. This would impede the existence of a reduced 
phase space. 

4. There is good evidence that the most general homogeneous cosmological model, the Mixmaster (or 


Bianchi IX) universe, is chaotic 67 69 72 


5. 


The well-known Belinsky-Khalatnikov-Lifshitz conjecture 73 posits that a generic (inhomogeneous) 
cosmological solution to the Einstein equations undergoes chaotic oscillations on approach to the 
initial singularity. 


6. It follows from a classification of generalized symmetries of the vacuum Einstein equations that, for 
closed spatial manifolds, no Dirac observables exist which are spatial integrals over local functions 
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of the metric and its derivatives 38 74 . This suggests that the vacuum Einstein equations are 
non-integrable. 

It may thus be no surprise that no Dirac observables are known for full general relativity - simply 
because they almost certainly do not exist in general in the first place. One should not get fooled 
and draw erroneous conclusions about gravitational Dirac observables and relational dynamics from our 
understanding of general relativity through simple integrable systems. 

Two important comments are in place. Firstly, this clearly does not imply that no general relativistic 
systems with Dirac observables exist. As mentioned in the introduction, general relativity subject to 
special asymptotic symmetries or special matter contents gives rise to Dirac observables. Furthermore, 
the simplest cosmological models are certainly integrable. But these systems do not account for a generic 
situation in general relativity. Secondly, weak non-integrability also does not entail the absence of gauge 
invariant observables. 

There are two notions in which gauge invariant observables appear: Firstly as non-differentiable or 


even non-continuous functions on the constraint surface (see also [24 25 for a brief discussion). This 
will be illustrated by means of toy models in sections [3] and [6] below. Thus, there is classically no 
conceptual problem, gauge invariant quantities exist by means of which one may interpret the dynamics. 
Yet, such generalized non-differentiable observables are not Dirac observables as they cannot be inserted 
into Poisson brackets and will thus fail to satisfy Poisson algebraic relations. 

Secondly, it is always possible (around a regular phase space point on the constraint hypersurface) to 
find sufficiently many smooth gauge invariant functions locally in phase space. That is, these functions 
will in general fail to be globally well defined or globally invariant. Again this will be illustrated in 
the examples in the main text. Indeed, the generalized Darboux theorem for constrained systems (see 
for instance 75 ) guarantees that locally there always exist sufficiently many phase space functions that 
commute with the constraints. 

Thus the issue we are discussing here is a highly non -perturbative one and results from global features 
of the phase space and constraint hypersurface. The drastic effect on quantization should not be a surprise 
as quantization depends strongly on global features 76 . This is again illustrated by the example, e.g. in 


section [3. 2 [ which differs from a fully integrable system by a periodic identification in the configuration 
space. 

The reader might wonder why we put much emphasis on general relativity and why weak non- 
integrability would not be an issue, e.g., for Yang-Mills theories. Certainly, depending on its constraint 
structure, a general constrained system will feature the same non-integrability issues. What is special 
about the canonical constraints in Yang-Mills-theories and the diffeomorphism (and even Gauss) con¬ 
straints in general relativity is that they are linear in the momentum variables (and their derivatives). 
However, linear constraints usually generate much simpler and integrable dynamics. The trouble maker 
for general relativity is the Hamiltonian constraint which is quadratic in the momentum variables and 
generates a much more complicated dynamics. Whence our focus on general relativistic systems. 

Lastly, we stress an essential conceptual and technical difference between unconstrained non-integrable 
systems and totally constrained, weakly non-integrable systems. The absence of differentiable constants 
of motion in the former does not pose a conceptual challenge: While constants of motion are convenient 
for fully solving the equations of motion, they are not strictly required for the interpretation of the system 
because one does not need to solve the dynamics in order to access the physical degrees of freedom. This 
also holds in the quantum theory where one does not need so solve the Schrodinger equation in order to 
access the physical observables; all observables on the Hilbert space are, in principle, ‘physical’. Whether 
the Schrodinger equation is difficult to solve or not is an unrelated problem. 

By contrast, in a totally constrained system one firstly needs to solve the dynamics in order to access 
the physical degrees of freedom - both in the classical and quantum theory. The ‘constants of motion’ 
of the constraints represent the physical degrees of freedom which are necessary for the interpretation 
of the system. Thus, the absence of Dirac observables has different repercussions than non-integrability 
in the unconstrained case. While classically this does not necessarily constitute a conceptual problem 
because the generalized observables are available, it potentially becomes both a conceptual and technical 









problem in the quantum theory since, depending on one’s method of quantization, one may not be able 
to construct a quantum representation of these non-differentiable quantities. This shall be one of the 
focuses of the main body of this article. 

3 Ergodic free particle dynamics 

In the remainder, we shall attempt to make some of this discussion concrete by studying simple toy 
models, rather than considering general relativity. Some of these models fail to be weakly integrable 
by being irreducible, but none of them will be fully weakly non-integrable. The advantage of working 
with models that are not fully weakly non-integrable is that they can be handled analytically in contrast 
to most chaotic systems, especially in a gravitational context (e.g., see the conceptual and numerical 
discussion in 3l]). Yet, by being irreducible, they already illustrate some of the behaviour generic to 
weakly non-integrable models. 

3.1 Free particle in the plane 

We begin with a trivial example: consider the Lagrangian with constant E 

L=^ x (x\+xl) + \E (3.1) 

as a function of ( X\,X 2 ) £ R 2 and A as well as their derivatives in ‘time’ t. This leads to the momenta 

dL m Xi dL 

Pi = = —r—, 7T A = -y=0, 

dXi A d\ 

where the second equation is a primary constraint. Using the Legendre transformations, the Lagrangian 
takes the form 


£ = i I p,+i 2 K -A||L + i|-£ 


and immediately leads to a secondary constraint (equation of motion for A) 


C = 


Pi 


pI 


2m 2m 


-B«0, 


(3.2) 


(3.3) 


where « denotes a weak equality. That is, A can be an arbitrary function of t. Clearly, C and n\ commute 
such that both are first class. Since the dynamics of A is trivial, we disregard it as a dynamical variable 
and simply treat it as a Lagrange multiplier - in analogy to the lapse in general relativity. This system 
thus has the interpretation of a free particle in R 2 with fixed energy. Given the arbitrariness of A as a 


function of f, it is clear that (3.3) generates reparametrizations in t such that t is a gauge parameter and 
not a proper time variable. 

The system is clearly both strongly and weakly integrable and p\ (or, equivalently, P 2 ) and L 3 = 
X 1 P 2 — X 2 P 1 constitute a set of independent Dirac observables parametrizing the reduced phase space. 
Choosing x\ as the relational ‘clock’ and denoting the values it takes by r £ R yields the relational 
observable 


, v P 2 , , . , x \j2mE-p\ L 3 

X2(t) = -(t-x i)+X 2 = sgn (p 2 ) -r- 

Pi Pi Pi 


(3.4) 


which, evidently, is invariant under (3.5). We shall now study how the integrability is affected by com- 


pactifying the dynamics on a torus and adding a few modifications to the constraint. 
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3.2 Compactification and anisotropic masses 

Next, we compactify the dynamics on a torus by identifying Xi +1 ~ Xi, 
anisotropic masses such that the constraint becomes 


C = 



p\ 

2 to 2 


-E 


0 . 


1 , 2 and, moreover, introduce 


(3.5) 


This system can also be viewed as two free particles on a circle and follows from a Lagrangian analogous 
to (3.1). The constraint (3.5) again generates reparametrizations in t. 


The constraint surface corresponding to (3.51 is compact and has topology 

C ~ / x T 2 , 


(3.6) 


where I = [—y/2miE, +\/2m/E] is a closed interval parametrized by pi (it could analogously be parametrized 
by P 2 ). The solutions to the equations of motion generated by (3.5) read 


X2 (t) = 


Pi . , \ 

—t + x 10 

mod 1 = 

Pi , , 

- 1 + Xio — 

Pi 

mi J 


mi 

.mi 

P2 , , \ 

- 1 + X 20 

mod 1 = 

P 2 l L 

—t + X 20 — 

P2 

m 2 J 


m 2 

-m 2 


-t + X 20 : 
) J 


(3.7) 

(3.8) 


where t is the gauge parameter along the flow of C, and [_-J denotes the floor function; here its values 
rii £ Z count the ‘winding number’ in 07 . For later purpose we note that 


, , Pi 

tan <p = 7 — 
P2 


(3.9) 


defines the angle of the trajectory on the flat torus, where 7 = m 2 /mi is the mass ratio (see figure [T]). 



Figure 1: A trajectory on the torus T 2 with angle given by (3.9 




Clearly, as before p± (or P 2 ) is a valid Dirac observable, parametrizing I and being constant along 
the trajectories. Since p\ strongly commutes with C, the model is strongly integrable and not fully 
weakly non-integrable. But does there exist another Dirac observable, F, in analogy to L 3 above, which 
is independent of pi and parametrizes the trajectories on T 2 where the non-trivial part of dynamics takes 
place? F(pi;xi,x 2 ) must satisfy at least one of d\F ^ 0 or d 2 F ^ 0 and, ideally, 


{F, C} = 8iF + ^ d 2 F : 

m 1 m2 


0 . 


(3.10) 


We shall show that no continuous function F on C exists which is constant along all trajectories and 
independent of p\. In particular, F satisfying (3.10) cannot exist. The model thus fails to be weakly 


integrable. To this end, we must distinguish two qualitatively very different cases: 
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resonant tori: 7 ^ £ Q such that all trajectories on T 2 are closed and periodic 


50 


51 


64 


non-resonant tori: 7 ^- ^ Q and the torus is filled densely by ergodic, non-periodic trajectories 


64 


50 


51 


Both the resonant and non-resonant tori fill the constraint surface densely since both the rational and 
irrational numbers fill / densely. However, the rationals have a zero Lebesgue measure in I such that 
the resonant tori are of Lebesgue measure zero in C. By contrast, the non-resonant tori are of Lebesgue 
measure 1 in C. 

It is the non-resonant tori on C which destroy the integrability of the system and the existence of a 
continuous F, independent of pi and constant along all trajectories. Namely, suppose 7 ^ ^ Q. Then all 
corresponding trajectories fill the torus densely, i.e., every neighbourhood of every point on the torus is 
intersected by each trajectory. Hence, F, for this fixed value of p\ has to take all its possible values it 
may assume on the entire torus in every of its open neighbourhoods. But this violates the definition of 
continuity. In fact, this situation holds for every point on the torus such that F is nowhere continuous 
on the torus. This models is therefore not weakly integrable. 

In particular, F fails to be differentiable on C such that one cannot insert it into Poisson brackets 
(globally on C ). Accordingly, there does not exist a sufficient Poisson algebra of independent Dirac 
observables parametrizing the solutions. Consequently, a reduced phase space does not exist. This 
immediately rules out quantization by the reduction method since we do not have Poisson algebraic 
structures available which we could promote into a commutator algebra in the quantum theory. 

Even more severely, the space of solutions (which for integrable systems constitutes the reduced phase 
space) is neither Hausdorff nor a manifold. The space of solutions is a quotient space S = C/ ~, where ~ 
corresponds to the equivalence relation y ~ y' <t=> ot^iy) = y' for some t, where y,y' £ C and a l c denotes 
the flow generated by (3.5): 


(i) <5 is non-Hausdorff: The equivalence relation ~ defines a quotient map Q : C —> S. The open 
sets Us in S are those which have an open set Uc C C as pre-image under Q. Denote by r pi a 
trajectory in the torus labeled by pi- Let pi give rise to a non-resonant torus. Since every t Pi 
densely fills this torus, every open set in C which is intersected by a given r pi is also intersected by 
every other r pi for the same p\. Indeed, there exist distinct r pi 7 ^ r pi : for example, {x\,X 2 ,Pi,P 2 ) 
and (; x\ + ^,X 2 ,Pi,P 2 ) with 7 ^ ^ Q do not lie on the same trajectory since no winding number ?ii 
exists which maps these two points into one another. Hence, every open set Us C S containing t Pi 
also contains every other r pi such that solutions on a non-resonant torus cannot be distinguished 
by open neighbourhoods in S. 


(ii) <5 is not a manifold: If S was a topological manifold, there would exist a homeomorphism fu : 
Us —> V C R n open for all open Us C S. In fact, n > 1, by degree of freedom counting (we have 
S ~ I x T where I is a one-dimensional manifold and T denotes the set of trajectories on a torus 
for fixed p±). Denote by Uc C C the pre-image under Q of Us, i.e., Us = Q{Uc)- Since Q is 
continuous by construction, we then have that fu °Q '■ Uc — > is continuous iff fu is continuous. 
Notice that fu o Q is constant along all trajectories in Uc- However, it follows from the above 
arguments, that there exist neighbourhoods Uc C C (intersecting non-resonant tori) such that no 
continuous fjy o Q : Uc — > K n>1 exists which is constant along all trajectories in i/crl But then for 
Us = Q(Uc) no homeomorphism /)-. : U —> R n>1 exists such that S is not a manifold. 


Therefore, the Marsden-Weinstein symplectic reduction fails in this toy model. 

It should be emphasized that this model is even fairly benign: (1) it is ‘semi-integrable’ since p\ is 
a Dirac observable, and (2) despite featuring ergodicity in the configuration manifold, the dynamics is 
not chaotic (trajectories either always stay parallel or have a linear relative distance growth in the gauge 
parameter). 

8 This would require at least two independent continuous functions Uc —¥ R which are constant along trajectories in Uc- 
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For a fully non-integrable - or even chaotic - system (and this is the generic case!), the situation will 
be worse: no differentiable (global) Dirac observables whatsoever will arise, impeding the existence of a 
reduced phase space and presumably any manifold structure for the space of solutions. Indeed, there is 
both analytical and numerical evidence that the space of solutions of the chaotic closed FRW model with 


minimally coupled massive scalar field features a fractal structure 66 70 


Nevertheless, as a generalization of Dirac observables, we can still have gauge invariant ‘observables’: 
gauge invariance means Constance along the orbits generated by the first class constraints and does not 
require differentiability. As seen above, in a generic case, such ‘observables’ will be highly discontin¬ 
uous functions without (Poisson) algebraic relations. (See section 5.6 for the construction of such an 


observable.) But at the classical level no conceptual difficulty arises; physical information need not be 
continuous. 

Alternatively, instead of global discontinuous observables, one may construct transient observables 
29-31 where temporal locality means validity only for sections of a trajectory. For instance, we can 


ask for the correlation of the ‘evolving’ aq with the ‘clock’ aq, in analogy to (3.4), despite this relation 


not constituting a global smooth relational Dirac observable due to the above arguments. To this end, 


it will be convenient to express the winding numbers of the trajectory, n*(t) := [^ht + aqoj in (3.7 
3.8), as functions of aq(t), aq(l), rather than t. We have to distinguish again between the resonant anc. 
non-resonant case. 

For a resonant torus, rq cannot be parametrized uniquely by aq(i),aq(f) (and initial data aqo) for the 
entire evolution because of periodicity. However, we can parametrize the rq for every period of the orbit 
by aq (t), X 2 (t) and aqo, £ 20 - Let N 2 be the number of complete revolutions in X 2 before a trajectory closes 
in on itself. Accordingly, since any fixed aq = r will be passed N 2 times during every period, aq(aq = r) 
will have N 2 solutions. Thanks to periodicity, these will also globally be all possible solutions. 

The situation is different for non-resonant tori. In this case, given initial data £ 10 , 2 : 20 , x i(t)> x 2 it) 
uniquely specify the winding numbers for the entire evolution, since the trajectories do not pass the 
same point twicej^] Thus, we may write rq(aq(t), aq (t); aqo, 2 : 20 ) which is implicitly well-defined for pairs 
(aq(t),aq(t)) lying on a trajectory defined by aqo, aqo- A numerical evaluation, however, would require 
infinite precision because every trajectory comes arbitrarily close to every point on the torus. There are 
infinitely many solutions for aq(aq = r): since every trajectory fills the torus densely, it will also densely 
cross a loop of fixed aq = r. The same applies to any closed loop on the torus such that the Gribov 
problem of fixing gauge becomes untractable in this model; for every gauge fixing, represented by such a 
loop, there exist densely many solutions. 

In either casej^Jwe can then solve the first equation in (3.8) for t and insert it into the second equation, 
while replacing the floor functions by the corresponding rq, to obtain the implicit relation 


1 P2 

2 : 2 ( 2 ") = (t - aq + nq(r, aq(r);aq, aq)) + aq - n 2 (r, aq(r); aq, aq), 
7 Pi 


(3.11) 


where r £ [0,1) is the value of the ‘clock’ aq, aq( t) is the relational observable yielding the value of 
aq when aq reads r, and aq,2q are the corresponding phase space variables. This relation between the 

9 More precisely, given x \0, X20 , x\ (/:). aq(t), the two winding numbers n; are found as follows: find 1 6 I such that 

Pi . 

Xi[t) — Xio = - 1 — 77 , 1 5 

m 1 

,,N P2 . 

X2{t) - X20 = - t — n2 

m 2 

gives ni,n2 E Z. There is at most one solution. Namely, suppose there exist distinct t' E M and ti^ti^ E Z also satisfying 
the last equations. Then one can easily check that we must have 

pi n' x - Tii 

7 — —-• 

P2 n' 2 — 712 

However, for a non-resonant torus 7— ^ Q which cannot be satisfied by 71^,71^ E Z such that 71^,712 cannot exist. Clearly, 
the case of one/no solution corresponds to x\(t), X2(t) lying/not lying on the trajectory with initial data £10, #20- 
10 In the resonant case, for a given period. 
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‘evolving’ X 2 and the ‘clock’ x\ is globally well-defined - albeit implicitly and multi-valued. While it 


is impossible to globally solve (3.11) for x-iij) (there exist ‘densely many’ solutions in the non-resonant 
case), locally along the trajectory, one can solve (3.11) explicitly: for every admissible pair of winding 
numbers (ni, 712 ), he. on every ‘branch’ of the solution, we find a linear relation between x-i and x-| (see 
figure For instance, on the ‘initial’ branch of the evolution we have rq = 0 such that, in analogy to 
(& 


X2 (t) = - —(t-X 1 )+X 2 . 
7 Pi 


( 3 . 12 ) 


We emphasize the temporal locality of these explicit relations; e.g., for an ergodic orbit, this relation is 
only valid on the ‘initial’ branch as no pair (xi,x 2 ) will be passed twice in the evolution. Furthermore, 
the range of permissible values of x 2 (t) and r is not the full T 2 and depends on the branch (711,712) (see 
figure [ 2 ]). 


x 2 , 

1 -T 7 - 

x 2 / 

0 xi 1 x\ 


Figure 2 : Each branch of a trajectory, given by an admissible pair (711,712), features a linear relation between x 2 
and x\. The range of permissible values of Xi depends on the branch; in the figure x\ £ [0,xi) and x 2 £ [x 2 , 1). 


Using rii(f) = [—t + Xio\ , it is not difficult to determine all values of n 2 compatible with a fixed value 
of 771 - and vice versa. In fact, for a trajectory on a resonant torus one can in this manner determine 
the finitely many admissible pairs of ( ni,n 2 ) for one period which, in principle, permits to explicitly 
compute all N 2 distinct solutions to ( 3.11| ). By contrast, for a trajectory on a non-resonant torus there 
exist infinitely many admissible pairs of (711,712), each giving rise to a distinct relation between x 2 and 
Xi, thereby preventing a global solution to (3.11). Nevertheless, locally we have explicit relations and 
predictions. 


4 Standard quantization 

A reduced phase space and a Poisson algebra of Dirac observables which could be promoted to a com¬ 
mutator algebra do not exist such that quantization by the reduction method is outright impossible. 
However, remarkably, Dirac quantization is possible. This model thus constitutes an extreme example 
for showing that reduced and Dirac quantization are in general not equivalent. Nevertheless, as we shall 
see shortly, also this standard Dirac quantum theory is ill and generically devoid of a semiclassical limit. 
In section [5] below, we shall discuss an alternative quantization which remedies this illness. 

11 Qualitatively, this is a gener aliz ation to continuous dynamics of the branch dependent relational observables of the 
chaotic discrete time Baker map 21 . 
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4.1 Dirac quantization and physical Hilbert spaces 

We choose wave functions tj)(xi,x 2 ) in the kinematical Hilbert space "Hkin = L 2 {S l x S 1 ) with xi,x 2 £ 
[0,1)- The momenta act as 


Pill) = -ihditp. 


(4.1) 


The momentum eigenstates 


tpki,k 2 (xi,x 2 ) = exp(27r ikixi) exp(2n ik 2 x 2 ) 

with (fci, k 2 ) £ Z 2 are a dense orthonormal basis set in Hkinj diagonalizing the constraint 

2mi 2m 2 

with discrete spectrum 


(4.2) 



E\(k!,k 2 ) £Z 2 }. 


Solutions to the constraint C 'i/’phys 


0 are defined by (hi, k 2 ) £ Z 2 such that 


fc? H—fc 2 = e := 
7 


717-1-E 1 

2 ( 7 rS) 2 


(4.3) 


(4.4) 


Solutions exist only for special (non-generic) If £ R + for which 0 will be in the spectrum. 

It is also clear that if ipkifa is i n the constraint kernel then so are i/)k 1 ,-k 2 i Moreover, 

if '*/ ; fci,fe 2 is i n the constraint kernel then 7 / 74 ,fcJ, can only be in the constraint kernel if 


7 = 


U 2 __ U2 

^ 2 ^2 

£/2 _ ^ 2 ’ 


(4.5) 


where 00 > 7 > 0. It follows that we have to distinguish the case 7 ^ 
For 7 ^ Q, there are only four cases: 


from 7 £ 


1. dim(ker((7)) = 0, the generic e for which no pair of integers ki exist that solves (4.4) 


2. dim(ker(C)) = 1, for the special case e = 0, where both ki vanish; 

3. dim(ker((7)) = 2, for special values of e = k 2 or e = - k 2 , when one ki vanishes and the other equals 

±/c; 

4. dim(ker(C')) = 4, for the special values e = k\ + - k 2 - 


Higher dimensional kernels do not arise since (4.5) cannot be satisfied for 7 ^ Q and ki, k[ £ Z. 


By contrast, the case of rational 7 is an intricate Diophantine problem, which possesses isolated values 
e with many solutions, e.g. for 7 = \ one can solve e = 9 as ( k\,k 2 ) = (3,0) or as (fci, k 2 ) = (1,4). 

We illustrate the discussion for the equal mass case 7 = 1 , where the problem is related to a well 
understood problem in number theory: In how many distinct ways w(n) can a given natural number n be 
written as the sum of squares of two integers ki,k 2 l We have to distinguish two cases, since w[n) regards 
k\ and k 2 as unordered, while in the constraint kernel k\ and k 2 have to be considered as ordered to 
appropriately count the dimension. Hence, if there is no integer k such that n = 2k 2 then the dimension 
of the kernel is 2w(n), otherwise it is 2(w(n) — 1). A classical result by Jacobi 
a given integer 


77 states that w{n) for 


n = 2‘ 




Pr 


9? 


(4.6) 
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where the Pi are the prime factors of n with p t mod 4 = 3 and the g, are the prime factors of n with 
g^mod 4=1, vanishes if at least one of the bi is odd (which is the bulk of integers n), and, if all bi are 
even, the number is 

w(n)= 2(ci + 1) • • • (c m + 1). (4.7) 

We see that the number of states in the constraint kernel is not bounded from above: There are al¬ 
ways exceptional values of e, such as e = 5 d for d £ N, for which the constraint kernel has dimension 
4 (d ± 1). For example, for e = 25 = 5 2 one has 4(2 + 1) = 12 states in the constraint kernel labeled 
by momentum eigenvalues (fci, k%) = (0, —5), (0,5), (—5,0), (5,0), (—3,-4), (3,-4), (—3,4), (3,4), 
(-4,-3), (4,-3), (-4,3), (4,3); or for e = 71825 = 5 2 13 2 17 one has 4(2 + 1)(2 + 1)(1 + 1) = 72 states 
in the constraint kernel labeled by the momentum eigenvalues (±15, ±400), (±32, ±399), (±76, ±393), 
(±81, ±392), (±113, ±384), (±140,±375), (±175, ±360), (±183, ±356), (±216, ±33), (±228, ±329), 
(±252, ±311), (±265, ±300). In conclusion, large constraint kernels exist but are rather exotic. 

In either case, the constraint kernel - and thus the physical Hilbert space H p h ys - is always finite 
dime nsional. H p h ys is spanned by the momentum eigenstates satisfying (4.4) and, since the spectrum 
(4.3) is discrete, a proper subspace of the kinematical Hilbert space L 2 (S L x S' 1 ), thereby inheriting its 


inner product. 


4.2 Physical observables 

Being finite dimensional, the physical Hilbert spaces are isomorphic to spin Hilbert spaces with observables 
on them taking the form of spin observables. In particular, no canonically conjugate pairs of Dirac 
observables exist on 'Hphys■ For example, the largest possible physical Hilbert spaces for 7 ^ Q are 
isomorphic to a two-qubit Hilbert space H p h ys — C 4 . This is a very small Hilbert space with few 
observables; a linearly independent basis set of observables on C 4 is comprised of 15 binary propositions 


and the identity (e.g., see 78 for a formulation in terms of questions). It can already be anticipated to 
be impossible to equip quantum theories on such Hilbert spaces with relational observables admitting an 
interpretation similar to that in the classical theory. 

Clearly, the momenta pi are observables on H p hy S , but also take the form of spin observables. The 
momentum eigenstates in H p h ys correspond to resonant tori for 7 £ Q and to non-resonant tori for 7 ^ Q 
since the ratio of the momentum eigenvalues is always rational. We shall exploit this shortly. However, 
we emphasize that, by means of linear combinations, one can engineer to be either rational or 

irrational, where (•, -) p hys is the inner product on Hkin restricted to the set of solutions to the constraint p 2 ] 

What about observables independent of the momenta? Classically, we had either global discontinuous 
or transient observables taking the place of Dirac observables involving positions. Unfortunately, the 
global discontinuous observables cannot be represented on H p h ys due to the discontinuity, while local 
observables of the type (3.12) cannot be represented as self-adjoint operators on H p h y s since they are not 
even classically globally defined. Leaving the issue of self-adjointness aside, one may hope that one can 


nevertheless give the local observables meaning - at least in the semiclassical limit as in 15 16 29 -31 
However, we shall now argue that this Dirac quantization of the toy model generically does not admit a 
semiclassical limit. Accordingly, relational dynamics cannot be meaningfully formulated; in particular, a 
relation to the classical relational dynamics cannot be established. 

Qualitatively, we seem to face the following issue: if there are few classical Dirac observables which 
can be nicely represented in the quantum theory, we have few physical states which can be distinguished 


12 For instance, suppose 7 = 1 and e = 2 k 2 . The physical state 

i/j = a exp (27 r i k(x 1 + X2)) + /3 exp( 27 r i k(x\ — X2)), 
can be easily checked to yield a ratio 

(Pi)phys 1 

(P2>phy S M 2 - Idl 2 

which, depending on a, /3 E C, can fail to be in Q. 
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by these observables. We shall circumnavigate this issue with a ‘polymer’ type quantization which admits 
the representation of more observables in section [5] 


4.3 A quantum illness: Absence of semiclassical states 


The generic failure of a suitable semiclassical limit is an immediate consequence of the previous discussion. 
Ideally, a proper semiclassical state should (a) be peaked on canonically conjugate Dirac observables, and 
(b) saturate the uncertainties in these Dirac observables. The trouble is, of course, that we do not even 
have canonically conjugate Dirac observables such that we cannot obtain such a semiclassical state in the 
first place. 

Notwithstanding, a finite dimensional constraint kernel and thereby an absence of canonically conju¬ 
gate Dirac observables on Hphys does not necessarily inhibit the construction of approximately semiclas¬ 
sical coherent states. For example, in the 2D oscillator model with fixed energy discussed in 16 30 the 


physical Hilbert spaces are likewise finite dimensional. However, there is a qualitative difference between 


that model and the present one: in the 2D oscillator model of 16 30 one has a reduced phase space and 


sufficiently many Dirac observables, one of which is analogous to an angular momentum with dependence 
on both position and momenta. These Dirac observables are even conjugate on the kinematical Hilbert 
space of which the physical Hilbert space is a proper subspace. One can peak on these Dirac observables, 
using the techniques in 79 , to produce sufficiently semiclassical elliptic coherent states for sufficiently 
large energies 30 . By contrast, in the present model, we do not have any Dirac observables involving 


the positions and generically there simply are too few physical states in order to build semiclassical wave 
packets. For such spin type Hilbert spaces, one would usually constructs semiclassical states using spin 
coherent states which, however, are not appropriate in the present case. 

Nevertheless, let us investigate what can be achieved by constructing ‘wave packets’. In general, one 
could imagine two distinct notions of semiclassical wave packets: 

1 . local wave packets: these are wave packets that are highly peaked near a (small) segment of a 
classical orbit, but which are allowed to decohere after finite internal time. Such local wave packets 
would be appropriate for the local relational observables - if existent. 

2 . global coherent states: these are wave packets that are required to be peaked along an entire classical 
orbit, which means that these states are not allowed to loose coherence in internal time. 

As before, let us choose x\ as the relational clock. We shall consider the initial wave function at 
X\ = 0 13 Furthermore, 


as in 


30 


we factorize the classical constraint as C = C- C +, where C± — 


Pi ± tj2miE — 4 p 2 generates evolution backward/forward in Xi time. The quantum constraint (4.4) 
then implies that the wave function satisfies a Schrodinger equation: 


■ h dip £ 
i h —— = H ip = s 

OX 1 


2miE - p 2 ip, 

7 


(4.8) 


where the operator s with spectrum ±1 arises as a quantization of the sign ambiguity of the spectral 


square root. The classical direction (3.9) is promoted to an operator 


tan (p = yp! p 2 1 . 

Clearly, the momentum eigenstates are eigenstates of the direction operator. The converse is not neces¬ 
sarily true: e.g., M + I3ip-k 1: -k 2 is an eigenstate of the direction operator but not the momenta. 
The physical interpretation of the tangent of the direction is the dimensionless speed of the effective 
evolution. 


We shall now use (4.8 1 to evolve initial wave functions into physical states resembling wave packets. 
To this end, we again have to distinguish the cases 7 € Q and 7 ^ Q. We begin with the latter. 


1 We are in the special case in which the space of initial wave packets is a tensor factor of 
functions is thus equipped with the inner product on this tensor factor. 


The space of initial wave 
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Figure 3: Example quasi-localized wave packet with k\ = 1, k 2 = 3. 
Left: Initial wave function at internal time x\ = 0. 

Right: Global evolution with the right-moving Hamiltonian. 


We exclude the case in which one or both ki vanish, since this implies that the initial wave packet is 
either completely delocalized or that there is no time evolution. We thus restrict to dim'Hphys = 4. The 
space of initial wave packets at X\ = 0 is two-dimensional with modes: exp(±27r* k 2 x 2 ). Thus, the most 
general initial wave function is of the form 


" 0 ( 0 , X 2 ) = a e 


2nik2X2 


■ 13 e 


— 2-Kik2X2 


Choosing s = ±sign(p 2 ) hr (4.81, H evolves the initial wave function to a solution of the constraint: 

0 ±( zi , 2 2 ) = ae 2ni{k2X2TklXl) + /3 e~ 2 ™ (k2X2TklXl ). 

It can be easily checked that independent of a, f3 


(Pl)phys _ fcl 
(P2)phys k 2 

The associated probability density reads 


(tan 0 ) p hys — +7 , 
K2 


(4.9) 


\tp±(xi,x 2 )\ 2 = \a\ 2 + |/3 | 2 + 23?[a*/3] cos ( ATrk 2 (x 2 + ) - 23[a*0] cos ( ATrk 2 {x 2 + — Xi) ) (4.10) 


It is straightforward to verify that the directions of constant |0 ±| 2 on the torus are (1,±||). That is, 
the evolution is coherent with dimensionless speed - independent of ot,f3. 

But this coherence is misleading. These states are not semiclassical at all: (1) the localization of 
these wave packets in position space is too weak to be interpreted semiclassically: The width of a peak 
in x 2 is 1 /fci and equals the separation of two subsequent peaks, so the ratio of width to separation is 
1. Good localization would require this ratio to be significantly smaller than 1 . This does not come as 
a surprise since there are no position Dirac observables. (2) It can also be checked that the spreads in 
the momentum Dirac observables are proportional to the energy which is not a semiclassical behaviour. 
Finally, (3) there is a mismatch in directions. We just saw that both the dimensionless speed and the 
ratio of the momentum expectation values are ±|+ while the expectation value of tan 0 differs by 7 ^ Q. 
This means that the peak of the ‘coherent wave packet’ does not actually follow a classical trajectory 
since the latter would have the direction determined by (tan0) p h ys - The wave packet is peaked on a 
periodic curve since ki £ Z, while a classical trajectory corresponding to these expectation values would 
be ergodic. 
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Figure 4: Example of an initially well localized wave packet with mass ratio 7 = 1 and e = 25. 
Left: Initial wave function at internal time x\ = 0. 

Right: Global evolution with the right-moving Hamiltonian. 



Figure 5: Example of a wave packet with mass ratio 7 = 1 and e = 125 that is initially localized on an 
orbit with periodicity 5. 

Left: Initial wave function at initial time x\ = 0. 

Right: Global evolution with the right-moving Hamiltonian. 


We illustrate these issues in figure [3] for a = which yields the states 


i>±(xi,x 2 ) = a cos ( 27rfc 2 ( x 2 =F fT x i 


Since we have exploited a full basis of Hphys, this is as good as it gets for 7 ^ Q. We thus conclude 
that in this case the quantum theory is sick and does not possess a semiclassical limit. We also note that 
this is the generic case because the irrational numbers are of Lebesgue measure 1 (in any finite interval) 
in the reals. 


7 € 


The discussion for values of e with small constraint kernels is completely analogous to the irrational 
case. Differences arise when the constraint kernel is larger. The simplest example with 7 = 1 and e = 25 
resp. e = 125, which is illustrated in figures [4] and [5| The figures show that it is possible to use the 
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allowed x 2 -frequencies 0,3,4, 5 resp. 2,5,10,11 to build initially quite well localized wave packets, which 
are either localized at a point (figure [4]) or at the intersection of an orbit (figure [5]). The price for the 
improved localization is decoherence: A coherent wave packet is of the form ip(x 2 A tan(^>)a;i), which 
would require a linear effective Hamiltonian w(fc) = ak in order to be generated from a general initial 
superposition 


*i>(o,x 2 ) = y, a i e2nik * X2 

i 


(4.11) 


where i runs over more than two values]^] (For superpositions with only maximally two of the a; ^ 0 
the same arguments as for 7 ^ Q apply.) Instead, the non-linear Hamiltonian (4.8 ) tha t appears in the 

solution of the constraint (4.41 yields ui(k) = 2miE— ~k? which evolves (4.11| into non-trivial 

superpositions 


ip(xi,x 2 ) = E /3j ipj(x 2 ± tan (</>j) Xi), 


with distinct <f>j which does not comprise a coherent wave packet on the torus p*] 

It follows from periodicity of the wave function and the nonlinearity of the effective Hamiltonian, 
that there can not be any well localized initial states that evolve coherently for long internal time. 
However, there are conditions for having a good initial localization and short term coherence. Good 
localization requires many allowed ^-frequencies that are small integer multiples of one another, while 
short term coherence requires that the Hamiltonian has a good linear approximation over the range of 
these frequencies. This requires that the allowed frequencies lie in an interval much smaller than \J1m\E. 
However, even such very special wave packets will loose coherence. The better the state is localized 
initially, the faster it will run apart. That is, global coherent states generically do not exist and local 
wave packets decohere very quickly. In fact, decoherence will generically happen on scales much smaller 


than the classical range of validity of a local relational observable of the type (3.12). As mentioned before, 


a quantum version of such local relational observables requires states which are coherent for sufficiently 
long. Thus, it is difficult to equip such local quantum observables with a tenable interpretation. 


4.4 Physical transition amplitudes and failure of semiclassicality 

Alternatively, one might try to probe the semiclassical sector of a theory by considering the transition 


amplitude between kinematical semiclassical sates 19 52 in order to evade the troubles of building 


semiclassical physical states. But, as we shall now argue, also such transition amplitudes do not admit 
semiclassical behaviour in the present case. 

Consider two wave functions |(a?i,pi)) and \{x 2l p 2 )), which are kinematical semiclassical wave func¬ 
tions peaked at (xi,Pi), then we consider 


IT(a;i,pi;f2,P2) 


_ ({x 2 ,p 2 )\P\{x 1 ,p 1 )) _ 

{{xi,pi)\P\{xi,pi)) ({x 2 ,p 2 )\P\{x 2 ,p 2 )) 


(4.12) 


where P denotes the projector onto the kernel of the Hamilton constraint. This is the timeless analogue 
of the time-dependent transition amplitude W(x±,pi, t\, x 2 ,p 2l t 2 ) = {(x 2 ,p 2 )\U(ti,t 2 )\(xi,pi)}, where 
U(t\,t 2 ) denotes the time evolution operator. 

The states are dynamically suitable semi-classical states with respect to the time evolution operator 
U if (1) the transition amplitude \W(xi,pi, t\\x 2 ,p 2 , t 2 )\ 2 « 1 whenever the time evolution of (xi,pi) 

14 A coherent wave packet is of the form P. L . af z e lnk P 2+olxl \ from which we read off that coherence requires u(k) = ak. 

15 By contrast, we emph asiz e that the coherence conditions for the relational evolution of the initial wave packets in 
the 2D oscillator model of 16| |30] is elliptic such that the non-linear square-root Hamiltonian of that model does produce 
approximately semiclassical states which do permit a temporally local relational dynamics [301. 
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from ti to f 2 is (x 2 ,P 2 ) and (2) \W(xi,pi,ti\ X 2 ,P 2 ,t 2 )\ 2 3 decays rapidljj^jas {x\,pi) moves away from 
the correct time evolution (x 2 , p 2 ). 

In analogy with this, one says that the states are dynamically suitable semi-classical states with 
respect to the projector P if (1) \W(xi,pi;X 2 ,P 2 )\ 2 ~ 1 whenever (xi,pi) and (cc 2 ,p 2 ) lie on the same 
classical orbit and (2) |W(aq,pi; 2 ; 2 ,jD 2 )| 2 decays rapidly when one moves away from an orbit. It is 
clear that the second requirement is insufficient for ergodic classical orbits, because an ergodic orbit that 
contains aq will come arbitrarily close to any af 2 and hence \W\ can not possess any position dependence. 
We will see this explicitly at the end of this subsection. 

The starting point for the semiclassical discussion is the fixing of a set of semiclassical sates. A partic¬ 
ularly versatile set of semiclassical states, which we will use here, are the so-called complexifier coherent 
states 96 98 . These are a straightforward generalization of the harmonic oscillator coherent states. The 


starting point for the construction of complexifier coherent states is the choice of an generalized Laplacian 
C, e.g. C = \pp. 

The operator C is called complexifier, because it is the quantization of a classical phase space function 
C, which is used to complexify the position variables x a using the definition 


E x 

j{x a , C} n , 
n! 


n=0 


(4.13) 


where {•, C} n denotes the ro-th nested Poisson bracket with C. E.g. for C = |p ■ p (with t the so called 
classicality parameter) one obtains the canonical complex phase space coordinates 

z(x a ) = x a +itp a . (4-14) 

The construction of the the complexifier coherent states is then straightforward: 

1 . One writes a formal wave function Sg 0 (x) that is completely localized at a point x 0 in configuration 

space as a formal sum 8g o = ^fe(*)l e fc) over a basis B. 

2. Define |(£ 0 )c) = e~ c ^ h S Xo := “0*0*0 e~ c ! h \ek). The sum \{x 0 )c) converges in the Hilbert 

space for suitable C. 

3. Consider \(x 0 )c) as a function of x a and replace x 0 with z(x a ). This defines the coherent state 


\(x 0 ,Po)c) ■= N\(x 0 )) 


for suitable C up to the normalization N. 




4. Normalize the states by setting N = {({x 0 ,p 0 )c\(x 0 ,p 0 )c)) 2 • 

For instance the usual harmonic oscillator coherent states can be constructed from S r 


= I dk 


n —ikx 


where |fc) denote momentum eigenstates, and the complexifier is given as C = ^P 2 with z(x) = x 
Following the recipe one obtains as complexifier coherent states 

\(x 0 ,Po)tp) = N J dke-i hk2 - ikx °+ kt r°\k) 

= Ne^ p ° f dke-^ {hk - p ° )2 - ikx 


3 \k), 

itp. 


1 k) 


(4.15) 


which are indeed the momentum space representation of the harmonic oscillator coherent states. 

We are now in the position to construct kinematic coherent states on a unit torus. We choose the 
standard complexifier C = |p 2 and obtain z(x) = x + itp. The construction then yields 

\{x 0 ,Po)c) = Ne^n p ° 2 | k), (4.16) 

kez 2 

16 E.g., the overlap function of two harmonic oscillator coherent states decays with the exponential of the phase space 
distance squared. 
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where | k) denotes a momentum eigenstate, now with discrete label k £ Z 2 . 

We can then derive the physical inner product between these semi-classical states as 

((fi,Pi)|P|(f 2 ,p 2 )) = M 2 Y. e -4^((Pi-2^fe) 2 +(P2-2rffe) 2 )-2«fc(ir 2 -a- 1 ) (4.17) 

fc£Z 2 \k\-\-k 2 ,/ 7=e 

where in A/" we absorbed the exponential factors eves A . (These factors will eventually drop out in W.) 
Here we used the projector onto the constraint kernel P = Sfcgz 2 -fe 2 +fe 2 / 7=e |^)(f*l 


withe = 2W' 


Any significant contribution to the sum will come from terms in which pj and p 2 are simultane¬ 
ously close to some 2tt hk a of the allowed k. If pj and p 2 are far apart, then there will be no k a and 
((£i,Pi)|-P|(:r 2 ,p 2 )) will be exponentially small. However, the two terms ((xi,Pi)\P\(xi,Pi)} may still be 
large, if each pj lies close to an allowed p. 

In the generic case all summands except for one will be exponentially suppressed because the allowed 
k do in general not lie close together. In this case the approximation gives 

e~ 4VR ( (p1 - 2n hk ° ) 2 + (pa - 27T hk„ ) 2 ) 

|W(fi,p' 1 ;f 2 ,p' 2 )| « e _^ ((gi _ 2fffifci)a+(g3 _ 2Tfifc2)a)) - (4-18) 

This will be exponentially small unless k 0 = k\ = fc 2 . In the later case |Wj ~ 1. This means that |Wj 
encodes information about whether pi and p 2 he close to the same allowed value 2nHk 0 or not. Moreover, 
we see that \W\ has no position dependence in the generic case. This lack of position dependence has 
a simple interpretation for ergodic orbits: The orbit [x\,p a ) comes arbitrarily close to any point x 2 . 
However, the present discussion applies to periodic orbits as well, as long as there is only one 2irhk 0 that 
is simultaneously close to p) and p 2 . Hence, there is a generic problem with the semiclassical behavior of 
W. 


In conclusion, the standard Dirac quantization of this model is ill because a semiclassical limit generically 
does not arise. We emphasize that this is qualitatively different from the relation between the classical 
and quantum theory in unconstrained chaotic systems where ergodic orbits can also not be traced out 
by semiclassical states 49 50 . Firstly, in many chaotic models one obtains so-called ‘quantum scars’ in 


the semiclassical limit which are increased probability densities along classically unstable periodic orbits, 
thereby at least approximately tracing out classical periodic orbits - in contrast to here. Secondly, the 
Hilbert space is usually infinite dimensional and separable which permits to build initially well-localized 
wave packets. Such wave packets will also run apart due to the sensitivity to initial conditions but one 
generally obtains at least a short time coherence which is more stable than in the present case. Thirdly, 
the conceptual consequences are different in this case: in the unconstrained chaotic case, there is no 
interpretational issue as one always has a unitary evolution with respect to an absolute time such that 
the fluctuations and decohering wave packets do not cause any difficulties in interpreting the dynamics. 
By contrast, in the present case, there is no absolute time and a relational interpretation of the dynamics 
is generically not possible with strongly decohering wave packets 29-31 . As seen, in the present case, 


even the more general interpretation in terms of transitions amplitudes is difficult. 

We emphasize that this model is still harmless: Part of our analysis was only possible because the 
momenta are Dirac observables. For instance, the direction operator is only meaningful thanks to these 
observables. The situation will be much worse for systems entirely deprived of global Dirac observables. 
The present analysis suggests that it might be impossible to meaningfully Dirac quantize such a system 
in this standard manner. Indeed, in section |6.1[ we shall see a model which only permits ergodic orbits 
classically and cannot be quantized in the standard Dirac manner, although not even being chaotic. In 
this case, without a quantum theory, the issue of absence of semiclassicality does not even pose itself. 


5 Introducing a discrete topology on the torus 

As discussed in section [2] and illustrated in the toy model in section |3.2[ the generalized observables are 
highly discontinuous in the case of an ergodic system. This suggests to consider an alternative topology 
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on the torus and with it an alternative measure. In particular, we consider a discrete topology on the 
configuration space (i.e. the torus). In such a discrete topology, every point is an open set such that the 
functions which are discontinuous with respect to the standard topology are actually continuous. The 
corresponding Hilbert space supports these functions such that we may hope to be able to represent the 
generalized observables on the ensuing (physical) Hilbert space - in contrast to the situation with the 
previous standard quantization attempts. 


5.1 Bohr compactification of the momentum space 


Such a discrete topology can be obtained by a so-called Bohr compactification of the dual space, i.e. 
the momentum space. Interestingly, the Bohr compactification appears as a quantization method in 


loop quantum cosmology 45 46 55 or polymer quantization 56 - 58 and also in the recent work on a 


new representation for loop quantum gravity 59 80 where the area operator resembles the constraint 
encountered here. 

The Bohr compactification leads to a representation of the canonical commutation relations that is 
unitarily inequivalent to the usual (Schrodinger like) representation. The Stone-von Neumann uniqueness 
theorem for quantum mechanics is evaded, as only the exponentiated momenta can be represented on 


the Hilbert space, whereas the momenta themselves do not exist as operators 81 


In the following we will describe the Hilbert space resulting from such a quantization procedure. We 
start with the position space representations, that is functions if{x\,x 2 ) on S 1 x S 1 . As before, the 
functions are periodic ip(x 1 + 1 , 22 ) = ‘f’{ x i,x 2 + 1) = if{x\,x 2 ). 

Equipping this space with a discrete topology means also that we will have an inner product based 
on a discrete measure 


(V’l 1^2) = dnd(xi,x 2 ) i>i{x 1 ,x 2 )i>2(xi,X2) 


(5.1) 


defined such that the states 


are orthonormal 


^ai,a 2 (£1,2:2) = s OL\ ,Xl 


(l^ai,Q2 ^a± 1 a , 1 ^Oi2 1 c 


(5.2) 


(5.3) 


Note that in (5.21 and (5.3) we use Kronecker delta functions (and not Dirac delta functions) defined as 

1 for a = a' mod 1 


firv.ry' — 


0 otherwise. 


(5.4) 


The states (5.2) form an (uncountable) orthonormal basis of the Hilbert space, which is given as the 
closure of the set of states of the form 


ip(x!,x 2 ) = S a i i>Xl S^ 2tX 


(5.5) 


where the sum is over finitely many indices i. Thus, the functions (5.5) would vanish almost everywhere 


(and would be highly discontinuous) with respect to the usual Lebesgue measure. 

Let us now pass to the momentum space representation, i.e. we consider functions if(k±, k 2 ) on Z x Z. 
In this momentum space representation the functions are known as almost periodic functions. The reason 


is that they are Fourier transforms of functions of the form (5.5). That is, we now consider (the norm 
closure of the set of) functions of the form 


tp(ki,k 2 ) = ^^c,; exp(— 2ni a\k\) exp(—2ni a‘ 2 k 2 ). 


(5.6) 
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The inner product for these functions is defined as 


I ^ 2 ) 


lim 

T —> 00 


1 

(2T + 1) 2 


^ ipi(ki,k 2 ) ip2(ki,k 2 ). 

|fci|,|fe 2 |<T 


and renders the (basis) states 

^ 01 , 02 (^ 1 ^ 2 ) = exp(— 2ni a\k\) exp(— 2m a 2 k 2 ) 


(5.7) 


(5.8) 


orthonormal. 

The transformations between position space and momentum space representation are given by 

4>{xi,x 2 ) = (V , xi,x 2 |V’}= lim .r.rr 1 , V' exp (2m X 1 ki ) exp (2m x 2 k 2 ) (/ fci , k 2 ) 

T—foo (21 + 1 “ z ' 

V ' |fci|,|fe 2 |<T 

ip(ki,k 2 ) = j dixd(xi,x 2 ) exp(—27r?'a;ifci) exp(— 2mx 2 k 2 ) ip(xi,x 2 ). (5.9) 

Note that functions of the form f(k±, k 2 ) = <5fc ljKl <5fc 2jK2 are n °t almost periodic functions and are not 
elements of the Hilbert space we are considering. In fact, these functions have zero norm in the inner 
product (5.7). 


5.2 Translation operators and their spectra 

This brings us to the question of how to represent the momentum operators - and hence the Hamiltonian 
constraint which is a function of the momenta. Momenta are usually quantized as derivative operators, 
however, we now deal with discontinuous functions. Accordingly, as mentioned in the beginning of this 
section, the momenta themselves do not exist as operators^] Instead, their exponentiations do exist and 
constitute well-defined translation operators which quantize the finite (symplectic) flow generated by the 
momenta. 

We therefore consider (in the position space representation) the translation operators 

(R%ip)(x 1 , 2 : 2 ) = i/>(xi+ n,x 2 ), (R%ip)(x 1} x 2 ) = i/j(x 1 ,x 2 + h). (5.10) 

This quantization thus introduces a novel (length) parameter /.i which warrants interpretation. The 
momentum operator p a , a = 1 , 2 is replaced by a /i-dependent operator 

PS = ( K ' 2 - R a M/2 ) (5.11) 

-2 

and the square of the momentum can be quantized as 

S a ■■= ^2 - 2I ) ■ ( 5 - 12 ) 

In order to keep the formulas in this section simple, we shall work with unit masses mi = m 2 = 1. 
However, there is no difficulty in working with non-unit masses and the results are qualitatively the 

17 This can be interpreted as a consequence of the uncertainty principle. The Hilbert space supports states that are 
sharply peaked in the position variables. Thus the uncertainty in the momentum would have to be unbounded, e.g the 
expectation values of (the square of) the momentum operator are not well defined with respect to all basis states. (The 
difference with the Dirac delta function state in the usual quantization scheme is that these are not elements of the Hilbert 
space and thus do not need to have finite expectation values with respect to observables defined on this Hilbert space.) This 
situation is resolved by the fact that the momentum operator is not defined on this Hilbert space. Only the exponentiated 
version exist, which results in a bounded operator and hence bounded expectation values. 
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same. For later use in section [6.1[ we mention at this stage already how distinct masses m a , a = 1,2 can 
be incorporated: 

:= 2^ - 2l) . (5.13) 

Alternatively, one may accommodate different masses (leading in general to a different spectrum) by 

{S'Ya ■■= ^- 21) • (5.14) 

2 p 2 m a 

The constraint we are considering (for this section with m\ = m 2 = 1) reads 

= Sf + S% - E. (5.15) 


Let us discuss the spectrum for this constraint and, more generally, for the (squared) momentum opera¬ 
tors. This will also elucidate the meaning of the parameter p. 

We start with discussing the properties - and in particular the spectrum - of the translation operator 
R%. This operator is invertible, isometric and therefore unitary. The spectrum for the translation 
operators has been discussed in detail in 59 in the context of constructing a new representation for the 


loop quantum gravity observable algebra, based on a BF-vacuum. There, the gauge group SU(2) is also 
equipped with a discrete topology such that one needs to replace the angular momentum operators by 
(group) translation operators. 

The spectrum of R% depends on whether p is rational or irrational. For this discussion, let us specialize 
to one copy of S' 1 and a translation operator R M ; the discussion directly carries over to the torus. 


p = p/q £ Q with p, q £ N and gcd (q,p) = 1: 


The Hilbert space is non-separable, however, the translation operator leaves finite dimensional sub¬ 
spaces invariant; e.g., a basis state ip a (x) = 5 a , x leads to the invariant subspace 


{(R^lpa = 'Ipa-Nfi I, N = 0 ,... ,q - 1}. 


(5.16) 


This makes it straightforward to diagonalize the translation operator on each subspace. The normalized 
eigenfunctions are given as 


v a , K (x) 


1 

7 ^ 


q -1 


z=o 


(5.17) 


with k £ Z. These have eigenvalues e 2 ’ rmM , thus k and k + q give the same value. (Also any a 1 = a + Ip 
with l £ 0,... ,q will give the same eigenfunction, but multiplied with a phase.) Hence, the spectrum of 
Rv reads 


{e 2 ^ | k = 0 , ...,q- 1}. 


(5.18) 


We can again start with a basis vector i[) a and define invariant subspaces, which now will be infinite 
dimensional (but separable): 

= clos(span{^) 0 ,_^ A( | N £ Z}). (5.19) 

(Remember that the labels are to be understood as modi.) Note that states ip a and ip a > with a = 
a' + Np mod 1 for some N £ Z define the same sub-Hilbert space. This defines an equivalence relation 
a' ~ a on the set [0,1) and we need to keep only labels in [0,1)/ ~. 
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The spectrum of R p on each of the is continuous and given by U{ 1). This is proven in more detail 
in 59]. The generalized, i.e. non-normalizable, eigenfunctions with generalized eigenvalue e 27r * p £ U( 1) 
and p £ [0,1) are given by 

U<*A X ) = Y^Sa+lKX ( 5 . 20 ) 

le z 


Replacing a with an a' ~ a from the same equivalence class leads to a multiplication of the generalized 
eigenvector with a phase. 

This result can be compared with the exponentiated momenta in the usual (smooth topology) quan¬ 
tization: There one obtains also a spectrum 17(1) = clos{exp(27ripfc)|fc £ Z} for irrational p. Yet this 
(seemingly continuous) spectrum is equipped with a discrete spectral measure dpa which arises from the 
discrete spectrum {k \ k £ Z} of the momentum operator. 

The compactification of the spectrum justifies the expression ‘Bohr compactification’. We can use the 


momentum representation (5.6) (adapted to one copy of Z) in order to display the generalized eigenvectors. 
These are given by (p is irrational) 


> 0 ) = 


^ g— 2nikoi^2ni l(p— /xfc) ^—2 / nikcx. 


' 6{p- pk) 


(5.21) 


where 5(p — pk) denotes the Dirac delta function with respect to a continuous measure dp on [0,1). Thus 
also here, the winding of the k £ Z values around the circle via exp(27 iikp) = exp(27r*p), leads to a 
continuous spectrum. 

Knowing the spectrum of the translation operator R p , we can construct the spectral decomposition of 
the Hilbert space %£. (Remember that this Hilbert space is separable, since the basis 
This decomposition is based on the identity 

/ dp e p u a ^p{x') 8 a -\-v ■ tAa-K'ju (■£)? 

J[0.1) 

where e~ 2mlp = (u atP \ ip a +i'p) and dp is the Lebesgue measure on S' 1 and coinciding with the spectral 
measure. Hence, for a vector ip £ R p we have the spectral decomposition 


5.19) is countable.) 


V = I dp Ua,p {u a , p | ip) 

J[ 0,1) 


into (generalized) eigenfunctions of the translation operator. This leads to a p-representation for states 
in HU as follows: 


V’ = c(l)ip a+ ip = dp F(p) 


(5.22) 


with 


F{p) = (u a , p \ip) ='$2 e 2 ™ lp c(l) and 

1 

c(l ) = {ip a +ip, | ip) = [ dp e 2nllp F(p) . 
“'[ 0 , 1 ) 

The inner product between two states can now be expressed as 

{ipi\ip2) = £ci(0c 2 (0 = [ dp F 1 {p)F 2 {p) . 

77Z J [o.D 


(5.23) 


(5.24) 
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Thus, we have constructed a unitary equivalence of the super-selection sector to L 2 (S 1 ,(dp)), in 
which the exponentiated momentum operators act as multiplication operators. The conjugated variable, 
i.e. x, is represented by the discrete variable l £ Z via x = (a + pi) mod 1. 

To summarize the different representations, we had started with a non-separable Hilbert space based 
on the discrete measure on the torus S' 1 x S' 1 . Via a Fourier transform we could change to a Hilbert space 
of almost periodic functions, depending on a momentum label k £ Z and with a (Bohr compactification) 
inner product (5.7). 

We then restricted to a separable Hilbert space on which the translation operator R} 1 has a 
non-degenerate spectrum. This allows to use the spectral decomposition of (for irrational p) to 
define the p representation, where p £ [0,1) parametrizes the eigenvalues exp(27rip) of the translation 
operator. The translation operator acts now as multiplication operator on an L 2 Hilbert space with 
compact configuration space S 1 and with the Lebesgue measure dp. 

Whereas the momentum variable has been (Bohr) compactified, the configuration variable x £ S' 1 is 
now parametrized by l £ Z, where x = (a + Z/z)modl. Note that due to the irrationality of p each l 
represents a different x £ S' 1 value, and hence there are no periodicity requirements in the variable l. 
Again we have a winding of ‘lattice points’ a + pi around S 1 , such that S' 1 is filled densely with these 
lattice points. 


5.3 Spectrum of and solutions to the constraint 

The Hamiltonian (or energy operator) for a free particle on a torus in the discrete topology therefore also 
depends on the choice of the parameter p and is of the form 


H* = (Sf + St) ■ 


The operator is a linear combination of unitary (and in particular bounded) operators. Therefore, 
the spectrum (and ‘kinematical’ energy) will be also bounded to be in the interval 


For p = p/q rational the spectrum reads 


2 p 2 


[ 0 , 8 ]. 


2 p 2 


{4 — 2cos(27tki^) 


2cos(27rK2/i)|Ki,K2 = 0,... ,q — 1}. 


(5.25) 


For q growing larger and larger the spectrum will fill out the allowed interval more and more densely. 
For p irrational we indeed obtain the full interval (5.25) as spectrum. 

The bound on the spectrum gives also an interpretation to the parameter p. It can be understood as 
the minimal step-size. Via the bound on the Hamiltonian it also determines the maximal energy to be 

Q fc 2 

given as • This is analogous to the situation in loop quantum cosmology, where the minimal step-size 


is related to the Planck scale and and the energy density is bounded by the Planck density 45 46 55 
Let us now discuss the solutions to the constraint 


C» = H» - E, 

where E obviously has to be inside the interval ^2 [0,8]. 
Solution to the constraint for p £ Q: 


The energy of H 11 applied to a state v ai ;Kl ® v a2tK , 2 is given by 


2p 2 


(4 — 2cOs(27TKi/i) 


2 C0s(27 TK 2 p )) . 
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Hence, for a generic E there will be again no solutions to the constraints . However, with growing q 
the allowed interval is filled more and more densely. Assuming a solution, labelled by (aq, a?, k±, K 2 ) 
does exist, there will be also a family of other solutions labelled by different a! x ,a! 2 - But these belong to 
different superselection sectors with respect to the translation operators R%. 


There is (at least) one degeneracy pattern, namely solutions with labels («i, K 2 ) lead also to solutions 
with labels (±ki,±K 2 ) (to be understood rnodg) and to (±/t 2 ,±«i) for |/ci| 7 ^ |« 21 - It could happen 
that for very special energies even more solutions arise, but we expect that these cases are even more 
non-generic than in the smooth topology case; finding solutions in the present case requires solving a 
‘trigonometric Diophantine problem’, instead of a quadratic Diophantine equation. 

As an aside, we note that allowing for states inside an energy interval E ± e rather than of fixed E 
in analogy to counting states in statistical physics - drastically increases the solutions for sufficiently 
large q. This is related to the notion of pseudo constraints and approximate (Gaussian) projectors rather 


than delta-function projectors 82 ■ 88 


The solutions to the constraints are normalizable, hence we can adopt the kinematical inner product 
as physical inner product. However, we face the same problem as with the quantization with smooth 
topology: Generically we have only four physical solutions. The difference is that the energy levels for 
which solutions do exist are now dense, but also bounded. 


Solution to the constraint for /.4 ^ Q: 

The spectrum of H M is continuous such that there will be (non-normalizable) solutions as long as the 
energy is in the allowed interval. H M has now generalized eigenvectors of the form 

£5 U a2 ,P 2 


with generalized eigenvalue 


2/i : 


777^ (4 - 2 cos( 27 rpi) - 2 cos( 27 rp 2 )) 


and p a £ [0,1). The (ai,a 2 ) label equivalence classes of angles, see the discussion after eq. (5.19). We 
can restrict to one choice of ( 011 , 0 : 2 ) - this means to restrict to one superselection sector with respect to 
the momentum operators and thus to a separable Hilbert space. 

The solutions for a fixed (allowed) E are labelled by several parameters: We can choose the energy 
0 < ei > E and as long as ei > E — 4 we can find a solution so that H^ = E. (We restrict to the 

generic case, i.e. exclude the boundary values 0 and §^2 for the energy E.) Additionally, there again 
exist sign degeneracies in the solutions. 

The crucial point is that we now have infinitely many solutions to the constraints. But as the solutions 
are non-normalizable with respect to the kinematical inner product, we have to construct a physical inner 
product in which these solutions become normalizable. 


5.4 Physical Hilbert space 

We now consider the case of irrational p and establish a physical inner product for the solutions of the 


constraint. Refined algebraic quantization (RAQ) or the group averaging procedure 89-92 is unfortu¬ 


nately not straightforward to apply as one would have to worry about a measure for the averaging of the 
constraints. Similarly, for the master constraint procedure 93 94 we would have to complete a direct 
integral decomposition of the Hilbert space based on the spectral measure for the constraint. 

We will also - as in the master constraint procedure - use the spectral measure, but define the inner 
product similarly to the RAQ procedure. 

To this end we use the p-representation discussed in section [A2j as the constraints are diagonal in 
this representation. That is we also restrict to one superselection sector (g> with respect to 
the translation operators Ri and R%. One can thus define a formal projector onto the solutions of the 


27 













constraints by 


P^ Q j = S (^E — cos(2npi) — cos(27rp 2 )^ (5.26) 

where E = 2 — 

This allows to define physical states as ‘projections’ of kinematical states F(pi,p 2 ): 

F p h ys (pi,P 2 ) ~ 5 (e - cos(27rpi) - cos(27rp 2 )) F(p 1 ,p 2 ) ■ (5.27) 

However as the ‘projector’ is given by a delta-function the resulting states are not normalizable in 
the (kinematical) inner product of T~L^ 1 ® . This is due to the fact that the projector property 

o P| a j = P{„} is n °t satisfied - one rather has an infinity appearing. 

To regulate this infinity one uses nevertheless the (formal) projector property and defines (in RAQ) 
the physical inner product as 


(FlI^Ws = (Fl I P{ a} M ■ F 2 } 

for the projections of two kinematical states]^] Here we introduced a measure factor A4(pi,p 2 ) with a 
positive function A4. 

We thus have for the physical inner product 


(Fi|F 2 )ph ys = J dpidp 2 M(pi,p 2 )5(2-- cos(2np 1 ) - cos(2np 2 ))F 1 (p 1 , p 2 )F 2 (p 1 , p 2 ). (5.28) 

We implement the delta-function in (5.28) by solving for p 2 

27rp^(pi) = arccos(i) — cos(27rpi)) 

where we defined pj to be the roots for arccos with 

p+e [0,1/2), P2 e [1/2,1), 

(that is p 2 {p\) = 1 — p 2 (pi)) so that we have sin(27rp + (pi)) > 0 and sin(27rp _ (pi)) < 0. We can then 
define 


f+(Pi) '■= F{p ll pt(pi)), f-{Pi) ■= F{p 2 , p 2 (pi)) (5.29) 

where for E £ [0, +2) we have pi £ [0, ^ arccos(£’ — 1)] U [1 — ^ arccos(£’ — 1), 1]. For E £ (—2,0) we 
have pi £ [ 7 ^ arccos(l — E), 1 — ^ arccos(l — E)]. (Here we take the arccos-root to be in [0, 7 r).) 

Note that the constraint equation 

cos(27rpi) + cos(27rp 2 ) = E (5.30) 

~ 2 

describes a deformed circle in (pi, p 2 ) space. E = 2—^E parametrizes the circumference of this deformed 
circle. This circumference is small for small energy E and vanishes for E = 0. This fact can be contrasted 
to the quantization scheme with continuous topology: There we had a constraint of the form hi + h 2 = E, 
with h\ and h 2 having discrete and positive spectrum. Thus one would in general expect that smaller 
energies lead to fewer states. This is indeed the case if one considers for instance two harmonic oscillators, 
where the physical Hilbert space is always finite dimensional and the dimension increases linearly with 
energy (as long as these correspond to energy values leading to any physical solution). In the discrete 
topology quantization of the particle on the torus this effect is now reflected in the phase space volume 

ls More precisely, we would need to specify a dense subspace in the kinematical Hilbert space on which the so called 
rigging map, i.e. the projector P/ ,. can be well defined as an operator from this dense subspace to a (algebraic) dual space 
of the kinematical Hilbert space. 
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of the constraint surface. However the physical Hilbert space is always infinite dimensional due to the 
spectrum of the Hamiltonians being continuous. In the standard quantization scheme for the particle on 
the torus the spectra are discrete, but the reduction of states is far more drastic than one would expect 
(e.g. from the harmonic oscillator example) and as discussed prevents a reasonable semi-classical limit. 

In summary, physical states are characterized by a pair of functions (/+,/_) (on the appropriate 
interval and with the correct periodicity behaviour in p). Having two such states (/+(pi),/_(pi)) and 
(< 7 +(pi), g-(pi)), the physical product is given as 


(9 I /) phys 


, M{p,p 2 {p)) 

J P 27r|sin(27rp 2 (p))| 


(g+(p)f+(p) + 9-{p)f-{p)) 


where we assume that M(p, p£(p)) = A4(p, p 2 {p))- 

As a test whether we obtain finite results for generic kinematical states, we compute the norm for the 

projections of the basis states V’q.ia := V’ai+Zip ( 8 ) V , ai-Hip- For the choice of measure factor M.(pi,p 2 ) = 1 

we obtain after a variable transformation 


1 2 I i’luh) = 


4 

47T 2 

4 

47T 2 


f f 1 

LL dxidX2 7^m^m 

f dx 2 1 

I^I- 1 


6(E — x\ — x 2 ) 


(5.31) 


This gives a finite result except for E = 0 (which one would exclude if one sticks with this physical inner 
product). The result simplifies if we choose M = 47r 2 | sin(27rpi)|| sin(27rp 2 )| to 


4 J J dxidx -2 5(E — x± — x 2 ) = 4(2 — |£j) , (5.32) 


which is finite for all allowed values of E. 


Remarkably, we have now obtained a physical Hilbert space in the form of an L 2 -space over a (mo¬ 
mentum) variable p. Note that there is the restriction for functions (/+,/_) to be of the form (5.29). 
Here the functions F from which the (/ + , /_) result are of the form (or arise as limit of functions of the 
form) 


f(pi, P2 ) = Y, c ^ l ^ e ~ 2nihple ~ 2nihP2 - 

h 


(5.33) 


Operators on the physical Hilbert space need, in particular, preserve this form of the functions 
(/+,/_). This is the case for the (exponentiated) momentum operators 

g27ripi ^ g 2Trip g2irip 2 ^ g 2TripJ(p) 

that act as multiplication operators in the p-representation. It also holds for the following derivative 
operator 

M (sin(2 7r p 2 )^--sin(27 r p i )^-) ~ sin(27rp^(p))^. 

(Here one uses the expression with p^{p) for f + (p) and p^{p) for /_(p).) This can be seen by evaluating 
the action of M on F(pi,p 2 ) explicitly: 


M ^ c(li,l 2 ) e ~ 2villPle ~ 2vil2P2 


2 


= VMluh + 1) - c(h,l 2 - 1)) - IMh + 1 , h) - c(h - 1, k))} e ~ 2nillpl e~ 2nil2P2 . 


Il,l2 
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Note that the operator M annihilates the constraint (which actually is of the form F(pi,p 2 ) in (5.33)) 

M (cos(27rpi) + cos(27rp 2 ) — E) = 0. 


Therefore M commutes with the constraint, which acts as multiplication operator in the p-representation. 
Hence, one can understand M as providing a configuration (Dirac) observable on the physical Hilbert 
space. It does not commute with the exponentiated momentum e 27Tlpi as 



— e 2mpi sin(27rp 2 ) 


which yields a non-trivial Dirac observable algebra. This allows us to construct eigenstates and semi- 
classical states that minimize the uncertainty relations with respect to these observables. 


5.5 Physical wave functions in the kinematical configuration space 

We have constructed a physical Hilbert space and Dirac observables acting on it. Let us discuss how the 
physical wave functions appear in the original parametrization. That is, we want to know how the states 
(P{Q}^ , )( a ' 1>a ' 2 ) l°°k likfi on configuration space. Note that the amplitudes of these wave functions do 
not have a priori a well-defined meaning (in the sense of not encoding directly a probability density for 
the position variables, as the latter are not Dirac observables). 

The question, of course, depends very much on the initial kinematical wave function ip that is projected 
to a physical one. Starting with a kinematical wave function 

V 1 = ^ ' C-{1\ , l 2 ) t/ , ai+ii/i,a2+i2A‘ 
h,h 

(where ip ai ,a 2 { x ij x 2 ) = S aiiXl S a2}X2 ), the projected wave function reads 

P{a}^ = ^ ] Cphys(h) ^ 2 )) 1pai+h^,a 2 +l 2 fi 

li,h 


with 

r n 2 

Cphy S {h,k) = dpxdp 2 M(pi,p 2 )8(2--cos(2irp 1 )-cos(2np 2 )) 

e 2 Tvi(h-I'^P! e 2iri(l 2 -l' 2 )p2 c n'i ^ 

:= 'EOih-i’uh-QW\A) 

h 5^2 

where 0 (Li,L 2 ) are the overlaps between the projections of kinematical basis states 

0(l\ - l\,h ~ I 2 ) = /dpid P2 .(2-^^- 003 ( 2 ^ 1 )-cos(27rp 2)) 

T. exp(27ri((/' 1 — h)pi + (1' 2 — h)p 2 ))- (5.34) 

^1 )^2 ’^2 

Figure [b] provides a plot of 0(Li,L 2 ) for E = 1 and with the choice M. = 47r 2 | sin(27rpi)|| sin(27rp 2 )|. 
Note that these overlaps 0(Li,L 2 ) are real due to the fact that the integral splits into a sum over four 
(sign sector) contributions, such that the imaginary parts cancel. 0(0,0) = 4 coincides with the physical 
norm of a (projected kinematical) basis state and also gives a significantly larger value than for the other 
(Li,L 2 ) ^ (0,0) values. 
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Figure 6: A plot of the overlap 0(L i, L 2 ). This also coincides with the plot of the physical wave functions 
obtained by projecting a basis state (peaked in the centre of the square) and not taking any winding into 
account. 

Let us consider the projection of a basis function with c(Z i,Z 2 ) = 3i 1 ,p 1 3i 2 ,p 2 - W.l.o.g. we can take 
(.Pi, P‘ 2 ) = (0,0). These are (maximally) peaked (kinematical) wave functions in position space. Thus, 
the uncertainty in momentum is also maximal (although technically not infinite, as the momentum is 
exponentiated to a unitary operator). To understand the plot of such a physical wave function, we observe 
that Zi, Z 2 give the position (aq,a: 2 ) of the value 0(Zi,Z 2 ) of the physical wave function where 

Xi(h) = (oq + hn) mod 1 , :r 2 (Z 2 ) = (a 2 + hp) mod 1. 

Accordingly, we need to wrap the values 0(Zi,Z 2 ) around the torus infinitely many times. Since p is 
irrational, no point will be visited twice such that no interference effect arises. This will lead in general 
to a highly irregular looking wave function in the original (kinematical) position space, e.g., see figure [7j 


This suggests that the projection of basis states will look in general irregular. A crucial point is 
here that the spread in momentum is maximal. The question arises if one can achieve more regular 
‘looking’ wave functions if one peaks in momentum instead. Here the difficulty is that the momentum 
eigenfunctions are generalized eigenvectors and thus not normalizable. As a result, we have to take into 
account some smearing. 

The translation operator eigenfunction (just for one copy of S l ) in position space are given by 

lei 

The set {(a + lp) mod 11 1 £ Z} will fill the interval [0,1) densely. To each point in this set one associates 
the amplitude a(Z) = e 2mlp . This can give again a highly irregular ‘looking’ function. However, one can 
choose p = up with «eZ. In this case we can consider 

e 2mKa u a ^) = 

zez 

which gives to nearby points in the set {( a+lp ) mod 11 1 € Z} ‘nearby amplitudes’. These cases substitute 
for the usual momentum eigenfunctions e 2nlK, t > in the standard quantization scheme. 
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1.0 



1.0 

Figure 7: A physical wave function with p = 0.051... taking 200 x 200 values 0{l\, Z 2 ) into account, hence 
around 10 windings. 


Switching again to the torus, we will have to satisfy the constraint 

E = cos(27rp!) + cos(27rp 2 )- 

Thus, for general E and a given p\ = n±p we can generically not expect to find a solution p 2 of the form 
p 2 = k 2 /t. Energies which do feature solutions (in that case at least four of those) with /Ci,« 2 € Z are 
given by cos(27t K\p) + cos(27 r/s 2 /i). Fixing one such energy, one will have solutions where pi and p 2 are 
of this form but also solutions where this is not the case. 

In consequence, ignoring non-normalizability, it is possible to peak on a solution of the form p\ = K\p 
and pi = K 2 P which will give a ‘regular’ looking wave function (assuming «i,k 2 are not too large). 
Unfortunately, wave functions peaked exactly at a pair (pi, p 2 ) are neither normalizable in the kinematical 
nor in the physical inner product (as the latter is a L 2 -Hilbert space with continuous measure in say 
pi). One should therefore take into account a smearing function E(pi,p 2 ) which needs to be of the form 
(5.33) (in particular periodic and smooth in (pi, P 2 ))- This will in general lead to an inclusion of p-values 
not of the form (ki/t, k 2 /x). 

It might still be possible to reverse engineer coefficients c(Zi,i 2 ) such that one obtains a ‘regular’ 
looking wave function. (In particular, it would be useful to know how (physical) eigenfunctions of L look 
like in position space.) But in any case these seem to be atypical. This makes an interpretation of the 
physical wave function in terms of the transient observables 29 - 3l], mentioned in sections 3.2 and |4j 
difficult and thus poses a challenge to relational dynamics. Instead, it seems advantageous to interpret 


the dynamics in terms of physical transition amplitudes 19 52 between kinematical states which are 


more generally applicable even when relational dynamics fails as in the present case. This we shall do in 
section 15.71 


5.6 Classical interpretation 

Before turning to the transition amplitudes let us provide a classical interpretation of the Dirac observable 


found in section 5.4 First let us comment that we use the term ‘(quantum) Dirac observable’, as the 


operator M commutes with the quantum constraint. As we will see, however, M corresponds to a highly 
discontinuous phase space function M (for which we did not use the term Dirac observable), that is 
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indeed constant along the gauge orbits. The ’trick’ that achieves the representation of M in the quantum 
theory is the de-periodification of the x variables into the l parameters. 

Let us first start with the classical dynamics. Having to work with exponentiated momenta, we 
replaced the constraint 

C° = \(jp\+pl)-E (5.35) 


with the ‘polymerized’ expression 


C* 1 = -=-(2 - cos(/xpi) - cos (pp 2 )) - E 


(5.36) 


approximates C° in the limit of small p (with the approximation improving for small momenta pi). 
However, as p cannot be represented on the discrete topology Hilbert space we rather quantize C^. There 
are therefore two different classical limits, one where one considers h to be small but p fixed and one 
where one additionally considers the limit of p —> 0. 


Let us consider p fixed. In this case, the classical trajectories with respect to the constraint (5.361 
are given by 


Xi{t) 

X2 (t) 


- sin(ppi)t + x w 

- sin (pp 2 )t + x 10 


mod 1, 
mod 1 


(5.37) 


All that changed with respect to the C° constraint is that the parameters pi are replaced by ^ sin (ppi). 
That is, the classical trajectories are the same if the momenta in the two systems are matched appropri¬ 
ately. Accordingly, the discussion in section [3)2] regarding the non-existence of smooth Dirac observables 
depending on the position variables still holds. 

Assuming for the moment that the position variables Xi are in R. instead of S 1 one would also replace 
the angular momentum Dirac observable L for the p = 0 case with 


pLv = sm(pp 2 )xi - sm(ppi)x 2 


(5.38) 


This is indeed a constant of motion for the non-compact configuration space, however with the periodic 
identification of the configuration space to the torus we rather have the time evolution 


pL^(t) = (sin(pp 2 )x w ~ sin(pp 1 )x 2 o) - {sm(pp 2 )m(t) - sin(/xpi)?r 2 (t)) 


(5.39) 


Whereas the first term on the RHS of (5.391 is constant in time, the second one is - via the winding 
numbers rii(t) - time dependent. Thus ZA is not a global Dirac observable (but a local or transient one). 
One thus has to modify ZA by adding the winding numbers to remove the time dependence 


M' := sm(pp 2 )(x 1 +n 1 (x,p)) ~sm(pp 1 )(x 2 + n 2 (x,p)) 


(5.40) 


Note that to define the winding numbers uniquely we need to select for each momentum pair (pi,P 2 ) 
and each orbit labelled by this momentum pair a reference point (ri,r 2 ) on this orbit that marks the 
position with vanishing winding numbers. ( ni(x,p ), n 2 (x,p)) are then defined as the number of windings 
along the trajectorjJ^] starting at (ri,r 2 ) and ending at [x\,x 2 ). (Changing the reference points leads to 
a momentum dependent (and therefore invariant) term that is added to M 1 .) 

M' is invariant along the trajectories as we managed to ‘unwind’ the configuration variables xi. This 
also applies to the parameter U, defined as x t = (ai + Up ) mod 1. If we translate the time dependence 

19 Note that the winding numbers would constitute an ideal clock in the sense of always going forward (here in steps) along 
the trajectory (although in steps of 1, which can be remedied by adding Xi to n,). However to ‘measure’ the corresponding 
clock one needs infinite precision in determining both positions Xi and momenta p,. 
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from Xi to li (and consider these variables for the moment to be continuous) the time evolution for these 
variables would satisfy 


h(t 2 ) - h(ti) _ sin(ppi) 


(5.41) 


h{ti) ~ h(ti) sin(pp 2 ) 

Thus the observable 

M := sm(np 2 )h - sin (ppi)l 2 (5-42) 

is invariant under time evolution. This is indeed the classical counterpart to the quantum observable M 


in (5.341. 


5.7 Semi-classical limit and transition amplitudes 

To probe the semi-classical limit (at fixed p) one can consider the so-called physical transition amplitudes 


19 52 


W(ip 11 ^ 2 ) := 


(V’l|V'2)phy S 


Vii 1 l|V’l)phys(t/'2|V , 2)phys 


(5.43) 


between kinematical states. That is, the physical transition amplitudes are given by the physical inner 
product (5.281 between (the projections of) kinematical states^] To probe semi-classical properties 
we can select as kinematical states semi-classical (or coherent) states, which are available more easily 
as compared to constructing proper physical and semi-classical states. A question one can consider is 
whether the transition amplitude between states peaked on phase space points lying on the same orbit 
is much larger than for states for which this is not the case. 

We work in the p-representation for the (super selection sector) Hilbert space <8> and thus 
with L 2 {S l x S' 1 , dpi d/ 92 ) ■ As semi-classical states we consider (a tensor product of) heat kernel or 
complexifier states on 17(1) x 17(1) 96-98 , which are defined a^ 1 ] 


fUp) ■■= E 


e 2 e 


n(EZ 


= • ( 5 - 44 ) 

V t \ 2 t t ) 

ne z x 7 

Here t is the classicality parameter and the states are peaked around po = a as well as lo = A/i (as a 
Fourier transform to the 1 picture shows). Furthermore, we applied in the second line of (5.44) Poisson 
resummationj^] This changes the sum in the first line, which for small t is slowly converging, to a fastly 
converging sum, where for sufficiently small t the n = 0 term may give already a good approximation. 
(This better convergence behaviour is however traded with a loss of manifest periodicity in p.) 

Let us now consider the physical inner product 


(■PerfjAJ ® ® F a i )}phys — 


dpidp 2 M(pi,p 2 )5 [e - cos(27rp!) - cos(27rp 2 )^ 

• (P 2 )- 


F t a°, K ^) Ft a °, X o(p 2 )F i 


cy\ ,A- 


(piK 


t 

a 2^2 ' 


(5.45) 


20 The literal interpretation as transition amplitudes applies in the situation where one can chose a clock variable T and 
if one chooses the kinematical states to include a factor %j)\ ~ 8{T — ri) and if )2 ~ S(T — T 2 ). This then gives the transition 
amplitude from clock time T 2 to n [95] . If such a clock variable T is not available the physical inner product gives a 
generalization of the concept of transition amplitudes; _ 

21 These are the same states as discussed in section [4 .4\ However, in |4.4| the U( 1) X U( 1) gives the (configuration space) 
torus, whereas here the U(l) X U( 1) is parametrized by the variables (pi,p 2 ) which are momentum labels. 

22 Poisson resummation uses (for suitable functions F) that F(ny/i) = ^ E] n ez w here F(k) = 

dy F(y) exp(— iky) is the Fourier transform of F(-). 
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We can proceed by solving the delta function for e.g. P 2 and are left with an integral over p\. If we use 


for the heat kernel states the Poisson resummed form in the second line of (5.44) we can apply a saddle 


point approximation for the integral, valid for the limit of small t. For this we consider only the n = 0 
term in each of the sums. 

It then turns out that the conditions for a saddle point (along the real integration contour) is that 
the labels ( a°, A°\(t 1 , A 1 ) describe two phase space points on the same orbit. That is, the n = 0 term is of 
the form 


2 ^ 
dpi -^-5 ± (pi) exp(-r ± (pi)) 


(5.46) 


where 


»(»■*) = 2((A?) 2 + (A£) 2 + (Al) 2 + (A*) 2 ) + 

- 4 ?r 2 [« - pi ) 2 + ( a° 2 - p 2 {pi)) 2 + K - Pi ) 2 + (02 - P2 (Pi)) 2 ] 


(5.47) 


and p 2 (pi) = ^ arccos (E — cos(27rpi)) denote the two roots for the arccos function. 

Assume that (cr°,<7 2 ) = (°'i> cr 2 ) =: (^li 0 ^) and that (a\,a 2 ) satisfies the constraint equation (for 
either the + or — branch of the arccos root). In this case we can easily identify the (global and local) 
maximum of ^(t^) to be given by pi = a\ and by choosing the same branch for p 2 as for a 2 . We can thus 
satisfy the saddle point condition for the real part of the exponent r. The condition on the imaginary 
part of r is, that it has to be stationary. This imaginary part is then given as 


3(7^) = 4tt[(A , 1 -A?)((7i-Pi) + (A^-A^)(c72-P2(Pi))] 


and its derivative is 


with 


dpi 


47T 


dp 


P 2 = 


—(Ai - A?) - (Aj - A ° 2 )—pi 

dpi 


sin(27rpi) 


sin(27rp 2 (pi)) 


We thus obtain the stationarity condition 


(Ai - A°) 
(Aj - K) 


sin(27rcri) 

sin(27T(j2) 


(5.48) 


(5.49) 


(5.50) 


(5.51) 


which indeed reproduces the classical trajectory (5.41). (Remember that the coherent states are peaked 
at l = A /t.) To be precise a saddle point approximation could possibly also be applied in cases where 
the cr and A parameters are not corresponding to a classical orbit. In this case, the maximum of the p 
and a dependent part of (displayed in the second line of (5.47)) will be negative whereas for the 

case that the momentum label are constant and on the constraint hypersurface, the maximum vanishes. 
Note that the A-dependent part of ^(r^) is canceled by the normalization of W. 

The expected peakedness property on the classical orbits can be confirmed by a numerical evaluation 
of the integrals. (We again choose M = 47r 2 | sin(27rpi)|| sin(27rp2)|.) We considered the approximation 
with (a) taking only the n = 0 term into account in all sums and (b) taking the n = 0 and n = 1 terms 
into account in all sums. (This refers to the sums after Poisson resummation.) For the choice of t = 0.01 
it turns out that there is no difference up to the order of 10~ 6 between the two approximations for all 
the transition amplitudes we computed. 

Figure [8] shows the norm of the transition amplitude 


W(a° = (0.2, p~2 (0.2)), A° = (; tLx(l° 2 ),tl° 2 ) \ a 1 = (0.2, p+(0.2)), A* = (0,0)) 


(5.52) 
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Figure 8: The norm of W as a function of with 1° determined by (5.53) and (Zj,Z|) = (0,0). 
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Figure 9: The norm of IT as a function of with 1° determined by (5.53) and (Zj,Zf) = (0,0). 


for E = 1 and t = 0.01. Furthermore the function L 2 

sin(27T(Ti) 


Li{1° 2 ) := 


sin(27rp2 (<7i)) 


(5.53) 


ensures that the coherent state is peaked along a classical orbit. (We can set w.l.o.g. l\ = l\ = 0.) By 
varying A 2 = tl% we display in Figure [8] the transition amplitude ‘along’ a classical orbit. 

As one can see the norm of W is (almost in this approximation) constant along the orbit. But the 
transition amplitude decays if we consider the transition amplitudes for states not satisfying the orbit 
condition ( 5.53[ ). Figure [ 9 ] shows the norm of W as a function of the deviation X defined by l° = Lifl^+X. 
We choose l% = 0 in the left and 1° = 30 in the right panel. In this sense we can probe the classical orbit 
with the physical transition amplitude involving (kinematical) coherent states. 


Note that once a choice of E and t is fixed (whose choice should reasonably depend on /i), the 
parameter /j, only appears implicitly. On the one hand in the translation of the l\ and Z 2 parameters 
into the actual positions X\ and X 2 and on the other hand in the identification of p with the momentum 
p via sin(27rp) = ^sin(/rp). Choosing p, sufficiently small, we can reconstruct the classical trajectories 
approximately for the (formally p = 0) original system with constraint C = ^{p\ + P 2 ) — E at least for 
a range of the momenta. Of course, we cannot compare to the standard quantization since as discussed 
the semi-classical limit for this case is unsatisfactory. 

Again, an essential mechanism that led to these satisfactory results is the ‘de-periodification’ of the 
X, variables into the li parameters. 
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6 Further toy models 


One can modify the original constraint (3.51 in a number of instructive manners. All systems in this 
section follow from a Lagrangian analogous to (3.1l and are accordingly reparametrization invariant. 


6.1 ‘Cosmological’ modification 

To render the constraint closer to the structure of a Hamiltonian constraint in relativistic cosmology we 
flip the relative sign between the momenta: 


2 2 

= _Pl_ _ JP2_ _ E 
2m\ 2m 2 


( 6 . 1 ) 


(We keep the dynamics compactified on a torus.) As before, the solutions to the equations of motion are 
given by ( |3.8| ), except that the second equation now features a minus sign in front of the p 2 . We have to 
distinguish (1) E = 0, and (2) E ^ 0. 

The E = 0 case is special as it implies 


7 Pi =pI 


( 6 . 2 ) 


and thus, as can be easily checked, if: (a) yq £ Q all classical solutions are periodic; (b) yq ^ Q all 
classical solutions are ergodic. 


(a) The argument in section 3.2 concerning the absence of a Dirac observable which is independent of p\ 
does not apply. Indeed, e.g. for 7 = 1, the following two functions 


Fi = sin ( 27 t(xi + x 2 )) (pi + P 2 ), 


F 2 = sin ( 27 t(xi - x 2 )) (pi - p 2 ) 


(6.3) 


are both smooth Dirac observables, independent of pi and commuting with (6.1) everywhere on C. 


Notice that Fi = 0 for pi/p 2 = —1, while F 2 = 0 for pi/p 2 = +1 such that both are required in order 
to parametrize the space of solutions, e.g. as F = T\ F 2 , which in this case will again take the form 
of a reduced phase space. Since everything is periodic, we can also write down relational observables 


as in (3.11) (replacing p 2 by —p 2 ), obtaining a finite number of solutions. This model clearly is weakly 
integrable (see section [ 2 ] for the definition of strong, weak and reduced (non-)integrability.) 

One can also quantize the model with the standard Dirac method. The solutions are given by 


yqfci = ±k 2 


(fci, k 2 ) £ Z~ 


(6.4) 


For yq £ Q there will always be solutions to this equation. In particular, for 7 = 1 one obtains an 


infinite dimensional separable H p h ys defined by (4.2) with 

(fci, k 2 ) = (0, 0), (fc, k ), (k, —k) with k £ Z\0. 


(6.5) 


This Hilbert space clearly possesses sufficiently many Dirac observables to build coherent semiclassical 
wave packets such that the issues of the quantization in section [4] do not arise. 


(b) Since all orbits are ergodic, the arguments of section |3.2| apply and we thus have no continuous 
gauge invariant observables with non-trivial dependence on the configuration space. The system is again 
strongly integrable but fails to be weakly integrable. 

Since there is no reduced phase space, a reduced quantization is again outright impossible. But now 


even a standard Dirac quantization is impossible as (6.4) will not have any solutions for yy ^ Q. 

However, a polymer quantization with discrete topology, as in section [5] will again produce a non¬ 
trivial quantum theory with solutions to the quantum constraint. Replacing p^/2m a by S% as in (5.13), 
solutions to the constraint 


C M = Sf - s% 
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are given by 
A 4 £ Q : 


« 1,«2 = 0 ,..., g — 1 , 

Pi, P2 £ [ 0 , 1 ) 


s.t. 

s.t. 


cos(2TTKifj,/y/mi) = cos(27t k 2 p / -\/ TO i), 
cos(27rpi/ v /m7) = cos(2n p 2 / y/rn^) ( 6 . 6 ) 


Clearly, there will always exist solutions to these equations. Alternatively, we can also replace p^/2m a 
by S'J'' in (5.14), resulting in the constraint 


(JlH _ _ g'f 1 


with a different spectrum and solutions defined by 

p £ Q : Ki, k 2 = 0, ..., q — 1, s.t. 7COs(27rKi/i) + 2(7 — 1) = cos(27tk 2 /x), 

p £ Q : pi,p 2 £ [0,1) s.t. 7COs(27rpi) + 2(7 — 1) = cos(27rp 2 ) (6.7) 

In particular, for p ^ Q we again obtain a continuous parameter set of solutions, resulting in an infinite 
dimensional non-separable 77 p hy S , which, however, will become separable once restricted to a superse¬ 
lection sector as in section [5] Qualitatively, this Bohr compactified quantum theory will be similar to 
the one in section [5] such that we shall not further discuss it in detail here. The striking conclusion is: 
this model cannot be quantized with the reduced or standard Dirac method, but does give an interesting 
quantum theory once resorting to a discrete topology. 


Finally, the E ^ 0 case is qualitatively indistinguishable from the model in section [372) Both resonant 
and non-resonant tori occur and the same arguments as before apply, showing that no Dirac observable 
which is independent of p\ exists as a consequence of ergodicity. Similarly, no reduced phase space 
exists. Instead, the space of solutions is non-Hausdorff and not a manifold. However, a standard Dirac 
quantization will be possible, albeit resulting in a troubled quantum theory as before. These illnesses can 
again be remedied using polymer quantization. 


6.2 Adding a free ‘clock’ 

Next, we can add a non-compact degree of freedom, acting as a ‘perfect clock’ (in analogy to the highly 


special free massive scalar field commonly employed in loop quantum cosmology [45 46 99)), such that 
only a subsystem will feature ergodicity. Consider the constraint 


c= _eL + j1 + j 1_ b> 

2 nit 2m\ 2 m 2 


E > 0, 


( 6 . 8 ) 


where t £ ffi. and 27 2 £ [0,1) as before. The constraint surface is of non-compact topology C ~ T> x K. x T 2 , 
where V is a disc of radius yf~E. The solutions to the equations of motion generated by (6.8) read 


f(s) 
xi(s) 
x 2 (s) 


Pt , , 

— s + to, 
m t 

Pi 

— s T X10 

777-1 

P2 . 

- S + X 20 - 

m 2 


Pi , 

— s + X10 I, 


L7711 

P‘2 

-m 2 


P 2 

— s + x 2 o I, 


(6.9) 


where s denotes the gauge parameter along the flow of (6.8) and |_-J is the floor function whose values 
n\,n 2 are the winding numbers in x\,x 2 as in (3.8). Clearly, we can choose pi,p 2 as Dirac observables, 


parametrizing V. Since these strongly commute with C, the system is strongly integrable. The non-trivial 
part of the dynamics thus takes place in the configuration space IxT 2 . We would like to be able to 
parametrize the trajectories by two independent Dirac observables. As in the model of section |3.2| we 
will have resonant tori for "/ p\/p 2 £ Q and non-resonant tori for r )P\jp 2 ^ Q. Consequently, in the latter 
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case, we will again have ergodicity in T 2 . However, the important difference to the previous case is that, 
thanks to the first equation in (6.9), we do not have ergodicity in K x T 2 . Consequently, the argument 


from section |3.2[ ruling out additional Dirac observables and a reduced phase space, does not apply to 
this model. Indeed, we can perform a global gauge fixing; e.g., the ‘clock’ T = p t t + piXi +p 2 x 2 is perfect 
in the sense that {T, C} ss 2 E > 0. Hence, T = const constitutes a global gauge which is intersected by 
every trajectory precisely once. Accordingly, we have a global Dirac bracket and thus a reduced phase 
space. Every smooth function on this reduced phase space is a global, Dirac observable. This model is 
thus also weakly integrable. 

Furthermore, it is trivial to write down relational observables with respect to either of the ‘clocks’ 
T, t. For instance, choosing t as the clock variable with values r £ R, we obtain from (6.9) 


xi( t) = 
x 2 (t) = 


(r-t)+x 1- — — (r-t) +xA, 
L nil vt J 


rrhp 1 ( 

mi p t 

m t p 2 , ,, , 

- (t — t)+x 2 — 

m 2 Pt 


m t pi | 

Lmi p t 

m tP2, ., , 

- (T ~ t) + X 2 

l m 2 Pt 


Thanks to the floor functions, these relational observables are discontinuous and thus not Dirac observ¬ 
ables. However, while non-differentiable in directions transversal to the constraint flow, these relational 
observables are differentiable along the flow, commuting with (6.8) 

A reduced quantization is, in principle, possible. Also, a standard Dirac quantization can be performed 
producing a quantum constraint with continuous spectrum 

t2 7.2 p.2 7.2 


r h 2 k,Q K 2 k 2 
l 2 mt 2 toi 


h z k‘ 2 
2 to 2 


- E 


k 0 £ K, (Aq,fc 2 ) £ z 2 j 


j^Z j^Z 

thanks to ko £ R. The kernel of the constraint will be parametrized by (/ci, fc 2 ) £ Z 2 such that ^7 + 


K*k. 


2^7 < E. Since the spectrum is continuous, one will have to construct a new physical inner product on 


the space of solutions. This can be done via a Klein-Gordon type inner product and taking care of the 
sign sectors (see, e.g., 100 ). Similarly, this model can be polymer quantized. However, we shall abstain 


from discussing this model in further detail. 


6.3 Parametrizing an ergodic system 

Lastly, we briefly consider a parametrized version of the two free particles on a circle, governed by the 
constraint 


C= Pt 


Pi 


JP 2_ 

2 to 2 


( 6 . 10 ) 


2rni 

where t £ R and Xi j2 £ [0,1) as before and C ~ M 3 x T 2 . Ergodicity in T 2 arises for ')Pi/p 2 £ Q- 
Nevertheless, this model corresponds to an integrable unconstrained system and thus clearly features a 
reduced phase space and Dirac observables. More precisely, t = t is a global gauge fixing, intersected by 
every trajectory. This defines a global Dirac bracket and thereby a reduced phase space. This model is 
thus even reducibly integrable. 

The relational observables with respect to the ‘clock ’t with values r, 


xi(r) = 
x 2 (t) = 


Pi 

mi 

P2 

m 2 


(t — t)+x 1 — 
(t — t) + x 2 - 


Py 

Lmi 

P2 

-m 2 


(r-t) +aqj, 

(t ~ t) + z 2 J , 


are discontinuous as before, yet constant along the flow of (6.10). 


Formally, a standard Dirac quantization of this model is equivalent to reduced phase space quantiza¬ 
tion, both resulting in a Schrodinger equation. The physical Hilbert space will be Hphys — L 2 {S X x S' 1 ) 
with sufficiently many Dirac observables. Clearly, this model can also be polymer quantized. 

23 Note that both the winding numbers ni,ri 2 are just constants for a fixed value of the parameter r as the argument of 
the floor function is constant. 


39 




















7 Conclusions and discussion 


We have pointed out technical and conceptual challenges that arise in the context of (weakly) non- 
integrable constrained systems for the concept of observables and, in consequence of that, for obtaining 
a healthy quantum theory which admits a reasonable dynamical interpretation. Weakly non-integrable 
systems are devoid of differentiable Dirac observables and thereby do not give rise to a reduced phase 
space. Clearly, gauge invariant quantities and the space of solutions as a topological quotient space always 
exist, however, in contrast to the integrable case, the gauge invariant observables are non-differentiable 
or even discontinuous, thus not admitting Poisson-algebraic structures, and the space of solutions will 
not carry a symplectic structure, nor will it in general even be a manifold. This makes a quantization 
by the reduction method outright impossible and poses challenges to a quantization through the Dirac 
and path integral method since the generalized discontinuous observables cannot be represented in the 
quantum theory when using a standard smooth topology. At the very least, this poses challenges to the 
interpretation of the quantum dynamics. 

This challenge is only seen if one takes the global features of the (constraint) system seriously. As 
mentioned, it is always possible to find sufficiently many gauge invariant functions locally around a given 
phase space point. It might be thus an effect which is not accessible by perturbation theory. 

Here, we have seen that the number of Dirac 


Global features influence quantization significantly 76 


observables is an indicator of the size of the physical Hilbert space and thus is an important litmus test 
for a reasonable semi-classical limit. Dirac observables and their properties are therefore important even 
if one can turn to alternative tools, like the physical transition amplitudes, in order to interpret the 
resulting theory. Another instance where global features turn out to influence heavily the properties of 
the quantum theory is the spectra of (Dirac) observables i°3- Here it is again the global flow of these 
observables (and thus also the global shape of the constraint hypersurface) which determines whether the 
spectra are discrete or continuous. 

In this article, we have considered this issue only for finite dimensional systems (e.g., for cosmological 
models) and illustrated them with a few simple toy models. The latter suggest that there is generically 
a serious problem in quantizing weakly non-integrable constrained systems using standard techniques 
and, in particular, a smooth topology. The ensuing quantum theory - if it exists at all - seems to be 
generically ill by not possessing a semiclassical limit that connects to the classical dynamics. We stress 
that this problem is qualitatively distinct from the typical technical challenges arising when quantizing 
unconstrained non-integrable or even chaotic systems 49-51 . In the latter case, one usually at least 


has a short term coherence of wave packets built from superpositions over different energies (although 
also there wave packets will quickly run apart due to the sensitivity of initial conditions) and classically 
periodic orbits survive as ‘quantum scars’ for eigenstates of sufficiently large energy. By contrast, our 
analysis suggests that weakly non-integrable constrained systems will generically not even have short term 
coherence, nor produce any classically periodic orbits (in the toy models regardless of the size of the fixed 
energy). In addition, also the physical transition amplitudes between kinematical states do not feature any 
semiclassical behaviour. These problems seem to be directly related to the scarcity of Dirac observables 
that can be represented in the quantum theory and, in consequence, to the resulting low-dimensionality 
of the physical Hilbert spaces. The key difference between the unconstrained and (totally) constrained 
case is that in the former one does not need to solve the dynamics in order to access the physical degrees 
of freedom, whereas in the latter case one does. The fact that these troubles are qualitatively distinct 
from the usual semiclassical subtleties of unconstrained chaotic systems is also manifested by the model 
of section |6.1| for which a quantization by standard methods is entirely impossible, despite the system 
possessing well-defined classical dynamics. 

We expect that the situation is qualitatively similar in the field theory case. Indeed, as discussed in 
section [2j there is good evidence that full general relativity constitutes a weakly non-integrable system 
such that generically differentiable Dirac observables and a reduced phase space may simply not exist rj 


24 The origin of the failure of integrability of the analysed toy models was ergodicity in configuration space. A well-defined 
notion of ergodicity usually requires compact dynamics which one would not expect for a typical gravitational scenario. 
However, non-integrability is more general than ergodicity or chaos and does not require compact dynamics. 
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To be sure, there exist integrable subsectors of general relativity (e.g., in the presence of special asymptotic 
symmetries or matter contents and in simple cosmological models), but these do not account for the 
generic situation. 

Classically, this does not represent a conceptual problem as one still has generalized gauge invariant 
observables at hand. Furthermore, as discussed in the main body, mathematically the relation between 
different dynamical quantities is still well-defined and gauge invariant even if the system is non-integrablc. 
Thus, one can, in principle, still interpret the classical dynamics relationally even if the relational observ¬ 
ables are (globally) only defined implicitly. In particular, one can resort to transient observables in the 
sense of [29||32] to have a local relational interpretation along sections of a classical trajectory. As seen 
in the toy model of section |3.2[ one can even explicitly solve for such transient relational observables on 
branches of the solution (see also 311 for a discussion of this in a chaotic cosmological model). Thus, 
locally predictions and explicit relations are possible and this better be the case as this is the situation in 
our own observable (spatial and temporal) vicinity of the universe. For example, also our solar system is 
chaotic 64 , but nevertheless for time scales relevant to us we can actually compute the orbits of celestial 
bodies to a high accuracy. 

However, in the quantum theory it is difficult to make sense of relational dynamics in the context of 
non-integrability - at least in the usual sense of locally scanning through the a priori timeless physical 
quantum state. For systems featuring the ‘global problem of time’, a given value of a relational ‘clock’ 
reading admits multiple classical solutions for the evolving degrees of freedom. States in the deep quantum 
regime will overlap on the classically distinct solutions and superpose internal time directions such that 
it is unclear how to interpret this relationally 


29-31 


This becomes especially challenging for systems 
where a given clock value admits densely many classical solutions. For example, in the discussed toy 
models the physical quantum states can be uniformly spread out (or, in the discrete topology, peaked 
almost everywhere) on the torus. It is not clear how to interpret this in terms of the relational paradigm. 
This is qualitatively similar to the analysis in 131 , showing that the effective relational dynamics in the 


closed FRW model with minimally coupled massive scalar field breaks down generically in the region of 
chaotic scattering. 

The example of a particle on a torus illustrates the interpretational challenge well: suppose that this 
model described a cosmological dynamics and that the time that the particle needs to traverse the torus 
once is a cosmologically large time, e.g. corresponds to one period in an (almost) cyclic universe. The 
problem of not having sufficiently many Dirac observables would only arise if one averages over the cycles, 
however not if one applies a more ‘local’ quantization approach. 

Do weak non-integrability and chaos thus constitute a fundamental obstruction for canonical quantum 
gravity? Not necessarily. While they pose challenges to quantizing general relativity and interpreting 
the resulting quantum theory, the present manuscript does offer a possible resolution. One may exploit 
the fact that the act of quantization is a rather ambiguous procedure. This, in particular, applies if the 
(physical) phase space features a non-standard topology, or as in the examples considered here is not 
a manifold or not even Hausdorff. It is important for gravity, as the physical phase space for general 
relativity will certainly not be given by a manifold. The key idea for the examples considered here 
is to modify the quantization procedure by refining the topology of the configuration space until the 
discontinuous generalized observables become continuous in the new topology and can be represented 
in the quantum theory. A similar strategy might work in the full theory (assuming an overwhelming 
amount of technical challenges can be overcome): one has to choose the quantization of the kinematical 
Hilbert space such that (a) the imposition of constraints leads to sufficiently many states allowing for 
a reasonable semi-classical limit and (b) (which as we showed is related to (a)), the resulting physical 
Hilbert space allows for a representation of (sufficiently many) gauge invariant variables. 

Here we have employed a discrete topology which leads to a Bohr compactification of momenta 
in the quantum theory. This corresponds to polymer quantization 56 - 58 which is the quantization 
technique of loop quantum cosmology 45||46{55 and of a new representation of loop quantum gravity 59 . 
Remarkably, polymer quantization has remedied the illnesses of the standard quantization techniques and 
produced a non-trivial quantum dynamics with sufficiently many observables in the toy models of this 
article. The price to pay for this is that the quantum theory becomes more challenging as one now 
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has to deal with highly discontinuous observables and states. But this is a technical problem and not 
a problem in principle. Furthermore, one has to deal with super-selection sectors and to interpret the 
additional quantization parameter which, however, in a gravitational context can be related to the Planck 
scale 45 46 55 . Note that this additional parameter arises from the necessity to use exponentiated 


momenta, which in turn can be understood as a consequence of the uncertainty relations and the choice 
to have states sharply peaked in the configuration variable as normalizable states. The latter choice is 
motivated by the desire to represent discontinuous observables in the quantum theory. 

While the present work provides non-trivial evidence that loop quantum gravity inspired quantization 
techniques give rise to an interesting quantum theory, it is of course an open question whether these 
techniques will also work for full general relativity - and ultimately whether our own universe is actually 
described by such a modified quantum theory. 

We emphasize that the challenges to a fundamentally relational interpretation of the quantum dynam¬ 
ics mentioned above remain the same for the polymer quantization. But this is not a dramatic observation 
because one may resort to an interpretation in terms of physical transition amplitudes between kinemat- 
ical states as proposed in 19 52 . Since the physical transition amplitudes are equivalent to the physical 


inner product, they are always applicable (whenever the physical Hilbert space exists) and indeed retain 
a meaning in the presence of chaos when the usual relational interpretation seems to fail. 

Before closing, we mention two alternatives to polymer quantization. Firstly, one may modify the 
constraint quantization technique by, instead of requiring the physical Hilbert space to be solely defined 
through the constraint kernel, also allowing states which are ‘almost’ annihilated by the constraints. One 
might attempt to motivate this operationally from the fact that it is fundamentally impossible to measure 
physical quantities exactly. This modification leads to the notion of pseudo-constraints and approximate 
rather than exact delta function projectors 82-88 . For instance, in the toy models the constraint has 


a form C = H — E. Rather than requiring that the energy be exactly E, one could permit values in an 
interval E ± e in analogy to counting states in statistical physics in small energy intervals rather than at 
exact energy values. Technically, this could be achieved by changing the constraint to C = e[H/e\ — E 
which generically leads to a larger kernel (and thus richer quantum theory) than the original C. While 
statistical properties of the constraint spectrum are unimportant when working with exact constraints 
(except whether the spectrum is discrete or continuous around zero), they will become useful for pseudo 
constraints. For example, in this case periodic orbit theory and concepts such as the Gutzwiller trace 
formula [50,102 may be applicable to determine the density of states in the ‘vicinity’ of the original 
constraint. 

Related to this, it is interesting to study a path integral quantization of the examples considered 
here. Path integral quantization for constraint systems is deeply related to Dirac quantization, as the 
path integral acts as a projector onto the solutions space of the constraints [3,[4]. However, this and the 
resulting solution space depends strongly on the choice of the measure in the path integral. For example, 
a standard procedure is to use a discretization of the path integral, which leads typically to a breaking 


of the (diffeomorphism) symmetry which is deeply intertwined with the (Hamiltonian) constraints 103 


The resulting path integral just imposes the pseudo-constraints discussed above, see also the discussion 
in 11041. More generally, one usually interprets the path integral as a procedure which requires less global 
information, and it would thus be interesting to see whether and how one can reproduce the results of the 
Dirac quantization discussed here. As the problem of having only discontinuous observables arises from 
averaging the (ergodic) dynamics over (infinitely large) time intervals, the issue discussed here might 
indeed be connected to ‘infrared issues’ that result from discussing minisuperspace path integrals and 
related probability interpretations for infinitely large universes (see also the recent works 105,106]). 

Lastly, it is worth mentioning that there exist modified gravity theories which do not face the problem 
of non-integrability as discu ssed in this article. For example, neither shape dynamics 107-109], 


nor 


Horava-Lifshitz gravity 110 111 , nor Einstein-aether theory 112 possess a Hamiltonian constraint p 5 ] 
Instead, they feature a proper Hamiltonian. Accordingly, in these theories one does not need to solve the 
dynamics in order to access the physical degrees of freedom such that chaos does not pose a fundamental 


25 See also 531 which suggest a version of uni-modular gravity to address this problem (at least for quantum cosmology). 
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challenge. 
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